Astronomy & Astrophysics manuscript no. HOD ' CFHTLS © ESO 20 1 2 

April 17, 2012 



Galaxy clustering in the CFHTLS-Wide: the changing relationship 
between galaxies and haloes since z ^ 1.2^ 

J. Coupon'-^'**, M. Kilbinger^^, H. J. McCracken^, O. Ilbert^, S. Arnouts'^, Y. Mellier^, 
U. Abbas^ S. de la Torre'^, Y. Goranova^ P. Hudelot^ J.-P. Kneib'', and O. Le Fevre^ 

' Astronomical Institute, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan 
~ Institute of Astronomy and Astrophysics, Academia Sinica, P.O. Box 23-141, Taipei 10617, Taiwan 
- - - ^ Excellence Cluster Universe, Technische Universitat Miinchen, Boltzmannstr 2, 85748 Garching, Germany 

CN ■ "• Universitats-Sternwarte Miinchen, Scheinerstr 1, 81679 Miinchen, Germany 

JT^ , ^ Institut d' Astrophysique de Paris, UMR7095 CNRS, Universite Pierre et Marie Curie, 98 bis Boulevard Arago, 75014 Paris, France 

>--^ . ' Laboratoire d' Astrophysique de Marseille (UMR 6110), CNRS-Universite de Provence, 38, rue Frederic Joliot-Curie, 13388 

Cn ■ Marseille Cedex 13, France 

^ ' ' Canada-France-Hawaii Telescope, 65-1238 Mamalahoa Highway, Kamuela, HI 9674 

CL( ** INAF-Osservatorio Astronomico di Torino, 10025 Pino Torinese, Italy 

■^^ ■ ' SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 

\^ • Received date / Accepted date 

ABSTRACT 

^^ ■ It has become increasingly apparent that studying how dark matter haloes are populated by galaxies can provide new insights into 

r \ ' galaxy formation and evolution. In this paper, we present a detailed investigation of the changing relationship between galaxies and 

^*^ , the dark matter haloes they inhabit from z ~ 1.2 to the present day. We do this by comparing precise galaxy clustering measurements 

(~ l . over 133 deg^ of the "Wide" component of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) with predictions of an 

^ f analytic halo occupation distribution (HOD) model where the number of galaxies in each halo depends only on the halo mass. Starting 

I ■ from a parent catalogue of ~ 3x10* galaxies at /^g < 22.5 we use accurate photometric redshifts calibrated using ~ 10* spectroscopic 

O , redshifts to create a series of type-selected volume-limited samples covering 0.2 < z < 1.2. Our principal result, based on clustering 

measurements in these samples, is a robust determination of the luminosity-to-halo mass ratio and its dependence on redshift and 

^ ■ galaxy type. For the full sample, this reaches a peak at low redshifts of Mj^'^'' = 4.5 x 10"/r'MQ and moves towards higher halo 






masses at higher redshifts. For redder galaxies the peak is at higher halo masses and does not evolve significantly over the entire 
redshift range of our survey. We also consider the evolution of bias, average halo mass and the fraction of satellites as a function 

^N^ ■ of redshift and luminosity. Our observed growth of a factor of ~ 2 in satellite fraction between z ~ I and z ~ is testament to the 

^ ■ limited role that galaxy merging plays in galaxy evolution for ~ IO'^/i^'Mq mass haloes at z < 1. Qualitatively, our observations are 

\j^ I consistent with a picture in which red galaxies in massive haloes have already accumulated most of their stellar mass by z ~ 1 and 

subsequently undergo little evolution until the present day. The observed movement of the peak location for the full galaxy population 

^O ■ is consistent with the bulk of star-formation activity migrating from higher mass haloes at high redshifts to lower mass haloes at lower 

f~^ ■ redshifts. 

f — ■ Key words. Cosmology: observations - large-scale structure of Universe - Galaxies: distances and redshifts - Galaxies: evolution - 

("^ ' Galaxies: haloes 

1 . Introduction tary network wh ere the most massiv e haloes lie at the nodes of 

• '^ • the cosmic web dSpringel et al.ll2006l) . Galaxies form as baryons 

V ■ In our current paradigm of galaxy formation, haloes of dark mat- f^jj j^j^ ^^e centres of dark matter haloes and the inflow o f cold 

5^ : ter grow from tiny imperfections in the early Universe. Against gas provides the fuel f or star formation (IWhite & Reeslll978l 

. C^ . the background (at the present day) of an accelerating Universe, gall & EfstafliiotJToSOh . While this scenario broadly reproduces 

structures grow "hierarchically": small haloes form first and j^e observed galaxy distribution over a large range of scales and 

merge to build up larger ones, resulting in a complex filamen- ^^^^^^ ^^^^^ ^^^^^^j ^p^^ questions remain concerning which 

processes regulate star formation inside haloes and how this 

Send offprint requests to: J. Coupon jg^^^ ^^ jj^g Hubble sequence observed at the present day. 

* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii .j^hese physical processes which drive galaxy formation and 
Telescope (CFHT) which is operated by the National Research Council i^^r ^- ^j» i^ii^^u 
,.,_„, ^„ A J . T ■ XT ■ ij o ■ J 1.TT ■ ri regulate star formation are expected to correlate closely to the 
(NRC) ot Canada, the Institut National des Sciences de 1 Univers ot the ° ,- , , , , , • , , , • rr^, ,- 
Centre National de la Recherche Scientifique (CNRS) of France, and "^^^s of dark matter haloes which host galaxies. Therefore, a 
the University of Hawaii. This work is based in part on data products Profitable avenue to pursue to better understand the physics 
produced at TERAPIX and the Canadian Astronomy Data Centre as of galaxy formation is to determine how the relationship be- 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- tween the stellar mass Mstai- and the dark matter halo mass M\, 
rative project of NRC and CNRS. changes over the lifetime of the Universe. In the local Universe, 
** e-mail: couponOasiaa. sinica. edu.tw abundance-matching studies have shown that at z < 0.1 the 
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stellar-to-halo mass ratio (SHMR, Mstar/Mh) reaches a maxi- 
mum of a few percent at Mh ~ lO '^Mp, much smaller than 
the universal baryon fraction ( iGuo et al. 2010). This indicates 
that the conversion efficiency of baryons to luminous galaxies is 
highest at these halo masses. In less massive haloes, supernovae 
winds are responsible for gas reheating, which reduces star for- 
mation and flattens the galaxy luminosity function at the faint 
end (Benson et al. 2003). In more massive haloes, feedback from 
active galactic nuclei dHopkins et al.l 120061 : ISomerville et al.l 
l2008l) can "quench" star formation, leading to formation of 
the "red sequence" of pa ssively evolving galaxies jBower et al.l 
120061: ICroton et alll2006l) . The relative importance of these dif- 
ferent processes as well as their dependence on environment, 
redshift and galaxy type is still very much an open question. 

Ideally we would like to make a detailed comparison be- 
tween the luminous content of dark matter haloes and the 
halo masses themselves as a function of redshift. Several 
methodologies exist to es timate halo mass es: for example, 
using satellite kinematics (iMore et al.l 20 III), X-ray tempera- 
tures ([Peterson & Fabian"2006': 'Cattaneo et al."2009^, or gravi- 
tational lensing (Heymans et al. 2006; Ma ndelbaum et al. 2003 
lAugeretafll20I0l: iLeauthaud et al.ll20I Ibl) . However, the range 
in redshifts over which these techniques can be applied is lim- 
ited, and amassing large samples is not always easy. 

At the price of making some assumptions concerning the 
profiles of dark matter haloes and their evolution in number den- 
sity, another way of inferring halo masses is from the observed 
galaxy clustering and mass functions. Usually, one makes the 
simplifying assumption that the number of galaxies in a given 
dark matter halo only depends on the halo mass, a n assumption 
which has b e en verified with A^-body simulations ( Berlind et al.l 
I2003L 120051: iBerUnd & Weinberd l2002t iKravtsov et alJ l2004 . 
Furthermore, the galaxy population of a halo, described by the 
"halo occupation distribution" (HOD), is assumed to be indepen- 
dent of environment and assembly history of the halo. Although 
some authors have shown that this "halo assembly bias" can 
have an observable eff'ect on galaxy clustering measurements, 
the magnitude of this effect is, f or the m oment, much smaller 
than the size of systematic errors (ICroton e t al. 2007). This will 
be an important effect to test in future surveys. 

In this paper we will use this model to investigate the chang- 
ing relationship of dark matter and luminous matter from z ~ 1.2 
to the present day. Until now, the majority of papers which 
have used HOD modelling to interpret galaxy clustering have 
analysed eith er large, low redshift surveys such as the SDSS 
(see, for e.g., IZehavi et al.ll20IOHBlake etani2008l:IZehavi et al.l 
I20IOI) . or smaller (~ I de g^) deep fields such as the COMBO- 
17 and COSMOS survey s jPhleps et al.l2006l:ISimon et al.l2009l: 
ILeauthaud et al.l l20IIbl) . For the moment, spectroscopic sur- 
veys at higher redshifts are also limited to small fields of 
view (IZheng et al.ll2007l:lAbbas et al.ll20I0h . lBrown et al.l ( l2008l) 
has interpreted clustering measurements in the framework of the 
halo model of slightly larger redshift baseline, although their sur- 
vey only covered ~ 7 deg^ and used five-band photometric red- 
shifts. None of these surveys had the required combination of 
depth and areal coverage to make reliable measurements from 
low redshift to z ~ I of large samples of galaxies selected by 
type and intrinsic luminosity. In contrast, this is an ideal task for 
the Canada-France-Hawaii Legacy Survey (CFHTLS) with its 
unique combination of depth, area, and image quality. 

The five-band u*,g', r', i',z' CFHTLS photometry allows us 
to measure photometric redshifts out to z ~ 1 .2 with ~ 3% accu- 
racy and small systematic errors ( ICoupon et al.ll2009l) . The large 
sky coverage, ~ 133 deg^ after masking, enables us to probe a 



large range of halo masses. In addition, from the four indepen- 
dent fields of the CFHTLS we can obtain a reliable estimate of 
cosmic variance from the data itself. In this paper we construct 
a series of volume-limited samples selected by type, luminosity, 
and redshift, and use the halo model combined with the HOD 
formalism to infer how each galaxy sample populates their host- 
ing dark matter haloes, and how this changes with look-back 
time. 

The paper is organized as follows. In Section |2] we present 
our data set and catalogue production, including data reduction 
and photometric redshift measurement. In Section |3] we out- 
line the techniques we use to estimate galaxy clustering, and 
in Section |4] we present our model. Results are described in 
Section|5]and discussed in Section|6] Conclusions are presented 
in Section|7] 

Throughout the paper we use a flat ACDM cosmology 
{Q,„ = 0.25, Oa = 0.75,//o =70kms-i Mpc"' andcrg = 0.8) 
with h - Hq/IQQ km s ' Mpc '. All magnitudes are given in 
the AB system. 



2. Observations, reductions and catalogue 
production 

2.1. The Canada-France-Hawaii Legacy Survey 

We use the "T0006" release of the CFHTLS Wid43 
dGoranovaet all 12009) . The CFHTLS Wide has the mean lim- 
iting AB magnitudes (measured as the 50% completeness for 
point sources) ~ 25.3,25.5,24.8,24.48,23.60 in u*,g',r',i',z', 
respectively. An important feature in T0006 compared to pre- 
vious releases is that for each MegaCam pointing ("tile"), the 
TERAPIX group provides a magnitude offset in each filter 
to account for tile-to-tile variations in magnitude zero points. 
These offsets are calculated using a stellar locus regression tech- 
nique (High et al,,i2009|). In each tile, objects are detected on a 
gri-x^ image dSzalav et al.l 119991) . Ga laxies are selected usin g 
SEXtractor "mag_auto" magnitudes dBertin & ArnoutsI 19961) . 
and colours are measured in 3" apertures which are all well 
matched to the median size of galaxies at our /' < 22.5 mag- 
nitude limit. The sky coverage of the four CFHTLS Wide fields 
is shown in Figure [1] 

The MegaCam /-band filter broke in 2006, and subsequent 
observations were made with a new /-band filter denoted by "y". 
Data from these two filters were treated separately with separate 
filter curves. Using a small data set observed with both the new 
y- and the old /-band, we detected no significant difference in 
photometric redshift accuracy. For the rest of the paper, we use 
/' to represent both filters. 

2.2. Photometric redshift estimation 

Our photome t ric re dshifts were measured using LePhare 
dArnouts et al. 20021) following the proc edures outlined in 
lllbert et al.l d2006l) and ICoupon et al] d2009l). Our prirn ary tem- 
plate set is composed of the four lColeman et aLJ dl980l) (CWW) 
observed spectra (Ell, Sbc, Scd, Irr) complemented with two 
observed sta r burst spectral energy distributions (SED) from 
iKinnev et al.l (1 1 9961) . This primary set of templates is optimized 
using the VVDS Deep spectroscopic sample. We performed an 
automatic calibration of the zero-points using spectroscopic red- 
shifts in the Wl, W3 and W4 fields. The calibration is obtained 



http://terapix.iap.fr/cplt/T0006-doc.pdf 



Coupon et al.: Galaxy clustering in the CFHTLS Wide 



-l 


' 1 




' 




1 ' 


■ 


' ' 








WDS deep 


■ 


-A 




















- 








1 










1 






















-6 








































-8 








































10 








































1? 


. 1 ... 1 ... 1 ... 1 .. . 


1 . 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- 


CFHTLS W2 


■ 
















■ 
















■ 












_ 














_ 


























: 


























- 








... 1 ... 1 ... 1 ... 1 . . 


- 



34 

Ra 



CFHTLS W3 
CFHTLS D3 
DEEP2 





























J^ 






^ 






^ ' 





6 


.11,11 


' 1 




' 1 




' 


■ 






WDS F22 




- 


4 


- 








- 














2 












































n 
















■ 
















■ 












2 


. . . 1 


. . . 1 . . . 1 




. 


- 



215 

Ra 



334 332 

Ra 



Fig. 1. Sky coverage of W1,W2,W3 and W4 fields. Each blue square represents one MegaCam pointing. Spectroscopic data (WDS 
F02, F22, and DEEP2) are in red. CFHTLS Deep fields overlapping with the wide survey (Dl and D3) are shown in green. 



by comparing the observed and modelled fluxes, and done itera- 
tively until the zero-points converge. For spectral types later than 
Sbc, we introduced a reddening E{B - V) = to 0.35 using the 
ICalzetti et aP (l2000h extinction law. 

Without near-infrared photometry, Lyman- and B aimer- 
breaks create degeneracies between high and low photometric 
redshift solutions. Following a Bay esian approach similar to 
iBenited (l2000l) and lllbert et all (l2006l) . we adopt a "prior" based 
on the observed spectroscopic redshift distribution to reduce the 
number of catastrophic failures. No redshift solution is allowed 
which will produce a galaxy brighter than Mg - -24. 

Our photometric redshifts were calibrated with spectro- 
scopic samples. The WDS Deep survey (Le Fevr e et al. 20051) 
is available in the Wl field. The WDS Deep is a pure 
magnitude-limited sample at /ab < 24 and we used 5 926 se- 
cure redshifts (flags 4 and 5) for photometric calibration, and 
for building the redshift distribution prior From the WDS F22 
("WDS WIDE", iGarilli et al. 2008) in W4 we used 4514 
galaxies limited to /ab < 22.5. This sample has a success rate of 
92%. Finally, in the W3 field, 1 267 redshifts from the DEEP2 



survey (iDavis e t al."2003') selected at /?ab < 24. 1 were used for 
photometric calibration. No spectroscopic redshifts were avail- 
able in W2. 

To estimate photometric redshifts, we adopted a slightly dif- 
ferent method than for previous studies to measure redshifts 
from the probability distribution function (PDF). Rather than di- 
rectly using the redshift that minimizes the x^ distribution func- 
tion, we chose the median of the PDF. It has been found that the 
maximum likelihood estimate can lead to a systematic concen- 
tration of photometric redshifts around discrete values ("redshift 
focusing"), although we checked that the effect is very small for 
redshift bins of size Az > 0. 1 . 

The redshift accuracy ("o"") of our samples with spectro- 
scopic redshift Zs and photometric redshift Zp is defined as 
caV(1 + Zs), where Az = |zp - Zsl represents the difference 
between spectroscopic and photometric redshift. We use the 
normalised median absolute d eviation as accuracy estimator 
(NMAD, lHoaglinetaDll983h . expressed as cta^ = 1.48 x 
median(|zp - zj). "Catastrophic" redshifts are defined as ob- 
jects with |zp - Zsl/(1 + Zs) > 0.15. The percentage of catas- 
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trophic redshifts is denoted by 77. The dispersion and the frac- 
tion of catastrophic redshifts increase from bright to faint sam- 
ples. The accuracy is o-a7(1 + zj = 0.034 (Wl) to 0.037 
(W4) at /' < 21.5 and reaches 0.039 (Wl) to 0.051 (W3) at 
21.5 < i' < 22.5. The failure rate increases from 77 ~ 1.4% at 
/' < 21.5 to 77 ~ 5.5% at 21.5 < /' < 22.5. Galaxies with redder 
colours (El, Sbc) have more accurate photometric redshift esti- 
mates (0-^7(1 + Zs) = 0.031, r] - 1.64%) than bluer ones (Scd, 
Irr and starbursts, o-^J{\ + z^) = 0.038, i] = 3.69%). 

We constructed a merged catalogue with unique objects from 
individual tile catalogues provided by TERAPIX. We identified 
duplicate sources positionally and chose the object closest to the 
centre of its original tile. Although photometric redshifts were 
estimated to i' < 24, we selected objects br ighter than i' - 22.5 
to maintain a low outlier rate (see Table 3 of ICoupon et al.l2009t) 
and to ensure catalogue completeness. This sample covers a total 
effective area of 133 deg^. We show the photometric redshift 
accuracy in the secure redshift range covered in this study, 0.2 < 
z< 1.2, compared to spectroscopic redshifts in Fig.|2] 

Stars appear as point-like sources on the image and can 
easily be identified using their half-light radius at bright mag- 
nitudes; however, at fainter magnitudes this estimator can be 
confused by unresolved galaxies. We refined our star-galaxy 
separation by combining this profile-magnitude criteria with a 
colour-based criteria. We computed the half-light radius limit 
for each tile from the median and dispersion of the brightest 
objects and classified as stars all objects brighter than /' = 20 
with flux radii limits smaller than this amount. In the magni- 
tude range 20 < /' < 22.5, we classified stars as objects smaller 
than the half-light radius and with;if^(star) < ;if^(galaxy), where 
;if^(object) is estimated from the stellar and galaxy template li- 
braries. 

We eliminated objects observed with less than three filters 
and with ;t'^(galaxy) > 100. We suspect that these objects are 
false detections or rare objects with wrong photometric redshift 
estimates. Most of these objects reside in specific areas, such as 
around bright stars and field borders where the masking proce- 
dure has failed. 

Using the photometric redshift, the associated best-fitting 
template and the observed apparent magnitude in a given band, 
we can directly measure the ^-correction and the absolute 
magnitude in any rest-frame band. Since at high redshifts the 
A:-correction depends strongly on the galaxy spectral energy dis- 
tribution, it is the main source of systematic error in determin- 
ing absolute magnitudes. To minimise ^-correction uncertain- 
ties, we derive the rest-frame luminosity at A using the object's 
apparent magnitude closer to /I X (1 H- z). We use either the m*, g' , 
r', /', or z' observed apparent magnitudes according to the red- 
shift of the galax y (the procedure is described in Appendix A of 
lllbert et aDl2005h . For this reason the bluest absolute magnitude 
estimate makes full advantage of the complete observed mag- 
nitude set. However, as the M-band flux has larger photometric 
errors, we decided to use M^-band magnitudes. 

3. Clustering measurements 

3.1. Sample selection and galaxy number density estimation 

Objects are selected with /' < 22.5 and divided into five redshift 
bins: 0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.8, 0.8 < z < 1.0, 
and 1.0 < z < 1.2. A large bin width (A, - 0.2) ensures a low 
bin-to-bin contamination from random errors {cr- < 0.1). We 
investigate effects from systematic errors in Sect. 13.21 The com- 
plete sample contains 2 924730 galaxies. To separate galaxies 
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Fig. 2. Photometric redshift accuracy in the 0.2 < photo-z < 1 .2 
interval. The figure shows the number density of photometric 
redshifts in Wl, W3 and W4, versus spectroscopic redshifts 
from WDS and DEEP2 surveys. The blue Une is a linear fit 
to photo-z versus spectro-z. 



into "red sequence" and "blue cloud" types, we used a criterion 
based on best-fitting galaxy templates. We define "red" galaxies 
as objects with best-fitting CWW templates (see Sect. 12.21 ) esti- 
mated as El and Sbc (early-type), and "blue" galaxies estimated 
as Sbc, Scd, Im, SBl and SB2 (late-type). Selecting galaxies by 
best-fitting template is a more robust way of selecting galaxies 
by type, as compared to a simple colour selection, and it is less 
sensitive to the effects of dust extinction. Our resulting colour 
selection in the absolute magnitude-redshift plane is shown in 
Fig. [3] We note that a small amount of red-classified galaxies 
at high redshift z > 0.8 have slightly bluer colours. Their dis- 
tribution peaks near the blue galaxy distribution, suggesting that 
these objects could be blue galaxies erroneously identified as red 
galaxies due to photometric redshift errors. We also note that 
a simple colour cut would not exactly reproduce our selection. 
However, in the interests of simplicity and clarity we keep the 
"red"/"blue" labels for the rest of the paper. 

We extract volume-limited luminosity-selected samples for 
each of the "full" (or "all galaxies"), "red" and "blue" samples, 
using Mg absolute magnitude thresholds (hereafter denoted as 
"luminosity threshold samples"), from Mg - 51og/! = -17.8 
(fainter threshold in the range 0.2 < z < 0.4) to Mg - 5 log h = 
-22.8 (brighter threshold in the range 1.0 < z < 1.2). The sam- 
ple selection is illustrated in Fig. H) For the rest of this paper 
we will refer to these samples as simply full, red and blue. Due 
to low numbers of pairs at small scales, luminous blue samples 
were discarded. We are left with 45 samples, each comprising on 
average ~ 153 000, ~ 70 000 and ~ 129 000 galaxies for a typi- 
cal full, red and blue sample, respectively. The sample properties 
are displayed in Tables IBTI IB. 2l and lB. 31 

Finally, in each redshift interval [Zmin;Zmax] we compute the 
galaxy number density: 



«g'l=Motal/ 



Q 






dz 



(1) 
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Fig. 4. Sample selection in the full (left), red (center) and blue (right) galaxy samples. In each panel, the galaxy number density in 
the plane Mgjz is shown; red rectangles represent the luminosity threshold samples. 



where D. represents the solid angle subtended by the survey, 
and dy/dz the volume element. Errors are estimated from the 
weighted galaxy number density field-to-field variance. 

3.2. Photometric redshift uncertainties 

Our modelled two-point correlation function is projected using 
the measured redshift distributions. In order to take into ac- 
count statistical errors on redshifts, we select galaxies in the 
redshift range considered and convolve the observed redshift 
distributions with the estimated photometric redshift errors, de- 
rived from the probability distribution functions. We construct a 
Gaussian error distribution for each galaxy centred on the me- 
dian redshift of the PDF, with a width corresponding to the 68% 
confidence limits of the PDF, and normalised to unity. We then 
sum these Gaussians to construct the redshift distribution resam- 
pled to a redshift bin width of 0.04. Redshift distributions for 
each sample are illustrated in Fig.|5] 

To further assess the quality of photometric r edshifts, we per- 
form t he cross-correlation analysis introduced in JBeniamin et al.l 
(I2OIOI) . The measurement of the angular correlation functions 
for galaxies in different photo-z bins is used to constrain the frac- 
tion of galaxies that are scattered into "wrong" redshift bins due 
to photo-z errors. We measure the bin-to-bin cross correlation 
function for the full, red and blue samples, respectively. A non- 
zero correlation between adjacent redshift bins is present in all 
cases. This may be due to the presence of large-scale structures 
extending over several redshift bins. More importantly, this is 
also due to photometric redshift scatter, which results in the leak- 
age of galaxies into neighbouring bins. This error contribution to 
the redshifts is taken into account by the convolution of the red- 
shift distribution with the errors, therefore, we do not consider 
adjacent bins further in the analysis of photometric redshifts un- 
certainties. 

The angular cross-correlation of galaxies in non-adjacent 
redshift bins is much lower, indicating a small fraction of catas- 
trophic outliers. We use the "global pairwise analysis" method to 
measure the contamination between two redshift bins / and j. In 
this approximation, the following linear combinations of the an- 
gular cross-correlation function Wij and the two auto-correlation 
functions, w„ and Wjj, respectively, are expected to cancel for all 
angular scales 9,, 



0. 



(2) 



d, = w, 



Ni 



N, 



m {fufjj + fufji) - W,(6,)j^fufjj - Wjj(er)^fjifu 



Here, A^, (A^^) is the observed number of galaxies in bin / (j). 
The contamination fij is the number of galaxies with true red- 
shift in bin ;, but misidentified into bin j, as a fraction of the 
true number of galaxies in bin /. For each bin pair (/, j), the leak- 
age of the other redshift bins is neglected. This a pproximation is 
valid for contamination fractions of up to 10% ( jBeniamin et al.l 
I2OIOI) . With this, the fraction of galaxies which stay in their bin 
is fa = 1 - fij. We fit the two parameters fij and fjj in Eq. |2] 
by performing ax~ null-test on d,. For the covariance {d,ds), we 
take into account the correlation between angular scales for each 
of the three correlation functions using a Jackknife estimate (see 
next section). We neglect the sub-dominant covariance between 
different correlation fun ctions. This correspond s to using the first 
three terms in Eq. A4 of lBeniamin et al.l (l20IOl) . 

Due to degeneracies between the parameters fjj and fji, large 
values for fij {i > j) cannot be ruled out in principle. However, 
for the full and red galaxy samples, the contamination fractions 
are consistent with zero in most cases. The blue galaxy samples 
are slightly worse, but contaminations are consistent with values 
between 2% and 10%. Together with the very low outlier rate for 
the sub-sample with spectroscopic redshifts (see Fig. |2]i, these 
results further strengthen our confidence in the photometric red- 
shift estimates in this work and to use them to measure angular 
correlation functions. Our measurements for the full sample are 
shown in Fig.|6] 

Bin-to-bin mixing caused by photometric redshift errors will 
reduce our measured clustering signal, given that objects in sep- 
arate bins are uncorrected. This would lead to an underestima- 
tion of halo masses, as low mass haloes tend to be less clustered 
than more massive ones. However, it is not trivial to predict how 
galaxy number density estimates would be affected by such er- 
rors, which may result in incomplete or contaminated samples. 
Since the halo number density decreases monotonically with 
halo mass, galaxy number density decreases with increasing host 
halo mass. Therefore, a contaminated sample will underestimate 
halo mass fitting estimates. Conversely, an incomplete sample 
will tend to overestimate halo masses. 

To understand how photometric redshift errors could affect 
our results as function of redshift and luminosity threshold, we 
compared our measurements in the CFHTLS Wide with those 
from the "Deep" component of the CFHTLS (iGoranova et al.l 
I2OO9I) . where m ore precise photome tric redshifts are available 
(see Table 3 in ICoupon et al.l 120091) . Photometric redshifts in 
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Fig. 5. Redshift distributions for the full (left), red (center) and blue (right) galaxy samples as function of absolute magnitude 
threshold. 



the CFHTLS Deep were computed using the method described 
in Sect. 12.21 and to avoid cosmic variance uncertainties we per- 
formed our tests in overlapping areas between the Wide and the 
Deep, Dl/Wl and D3/W3, over a total area of 2 deg^. 

These tests show that faint luminosity threshold samples are 
slightly incomplete at low redshift, probably as a consequence of 
catastrophic errors moving low redshift galaxies to higher red- 
shift. We report that the difference in galaxy density between 
the Wide and the Deep estimates rarely exceeds the field-to-field 
variance error estimate. However, brighter and brighter samples 
become more and more contaminated by spurious objects, and in 
the worst case (the brightest sample in the range 0.8 < z< 1.0), 
the number of objects in the Wide is more than three times higher 
than in the Deep, much larger than the 40% field-to-field vari- 
ance. A higher galaxy number density estimate will enhance the 
effect of a reduced clustering signal, which may result in under- 
estimated halo masses for brighter samples up to a few sigmas. 
In the worst case, this represents a bias of ~ 0.5 in log M^in. 

3.3. Angular correlation function 

We measure the two-point angular correlation function w using 
the lLandv & Szalavl(ll993h estimator 



w(e) 



NAN, - 1) {DP} M- - 1 (PR) 
NdNi-l){RR} " Na (RR) 



+ 1. 



(3) 



where (DD) is the number of galaxy pairs, (RR) the number of 
random pairs, (PR) the number of galaxy -random pairs, all in a 
bin around the angular separation 6. Nd and A^i are the number 
of galaxies and random objects, respectively. A random cata- 
logue is generated for each sample. In order to find an acceptable 



compromise between small Poisson errors and computational re- 
quirements we scale A^i depending on the number of galaxies in 
our input catalogues, from 50xA^d in the case of low density data 
catalogues to 2 x A^d for high density catalogues. We measure w 
using a fast two-dimensional tree code. At large separations 6, 
instead of counting individual galaxy pairs, we correlate boxes 
by multiplying the number of objects in the two boxes. For that, 
we define a threshold angle a^, - 6/db, where db is the box size. 
Below this threshold, the number of objects inside the box is 
taken into account instead of individual objects. We found that 
ff,,, = 0.05 gives accurate results at low computational cost. 

We correct our galaxy correlation meas urements for the "in- 
tegral constraint" (iGroth & PeeblesI Il977h . a correction factor 
which arises from the finite area of our survey Q^: 



W{0) = Wnies(6') + Wc ■ 



(4) 



The correction factor wc can be estimated as follows: 

wc^-^ ffw(0)dnidil2, (5) 

assuming a simple power law fitted on the data with slope y and 
amplitude A,,.: 

Wn,es(0) = A,,e'-^ -Wc^ A^O'-'' ' Q . (6) 

As A,t varies with each sample, we first compute C using Monte 
Carlo integration over random pairs 



C 



2 e^-yRR(ff) 



ZRRie) ' 

and then correct Wmes(0) using 

ffl-7 
W{e) = Wmes(0) 



e'-y-c 



(7) 



(8) 
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Fig. 3. Top: type selection based on best fitting templates for the 
"red" sample (left) and the "blue" sample (right) in the CFHTLS 
Wide. We show the colour distribution (M„ - Mr) as function 
of absolute magnitude (Mg) and redshift (top to bottom) and 
number counts for "red" and "blue" objects on the right panels. 
Bottom: number of red and blue galaxies as function of magni- 
tude and redshift. 



The integral constraint for the four Wide fields is C ~ 0.5. 
Assuming y - 1 .8, this leads to a correction of ~ 10% at = 0.1 
deg and up to a factor of two at = 1.0 deg. We note that as- 
suming a power-law for w{ff) to estimate C could be a source of 
error; it is neglected in the current analysis as the correction is 
smaller than our Jackknife error estimates. 

We combine our galaxy samples from the four CFHTLS 
fields into a single catalogue for the correlation function mea- 
surement. Computing w{ff) on the four fields independently, and 



Fig. 6. Cross correlation analysis between redshift bins for the 
full sample. Top: auto-correlation in the redshift bin 0.4 < z < 
0.6 (straight line) and cross-correlation between the bins 0.4 < 
z < 0.6 and 0.8 < z < 1.0 (dotted Hne) and between 0.4 < 
z < 0.6 and 1.0 < z < 1.2 (dashed line). Bottom: quantitative 
estimates of the contamination (percentage of galaxies scattered) 
from a pairwise analysis between redshift bins 0.4 < z < 0.6 and 
0.8 < z < 1.0 (left) and between 0.4 < z < 0.6 and 1.0 < z < 1.2 
(right). The contours show the 68.3 (blue), 95.5 (green) and 99.7 
(red) confidence regions. 



using a pair-weighted average would not lead to the same clus- 
tering estimate at large scales. Relative photometric offsets are 
expected to vary more from one field to another than within the 
same field, as the four fields are non-overlapping. Since our par- 
ent catalogue is selected by apparent magnitude cut (/' < 22.5), 
even a small difference in photometric offsets will lead to a spu- 
rious field-to-field variation of the galaxy number density. In 
addition, different photometric offsets could bias the photomet- 
ric redshift selection and also increase the field-to-field density 
fluctuations. A 0.02 magnitude offset variation would result in a 
~ 2% variation in galaxy number density, which may bias the 
two-point correlation at very large scales. Because we expect 
these effects to occur at the scale of individual field (a few de- 
grees), we adopt a conservative cut and limit our measurements 
to 6 < 1 .5 deg. We do not measure w at separations below 3'.'6 to 
avoid blended objects. Two-point correlation function measure- 
ments and errors are provided in the appendix ICl 



3.4. w{0) error estimates 

We estimate statistical errors on the two-point correlation func- 
tion using the Jackknife internal estimator. We divide all sam- 
ples into A^ = 68 sub-samples of about 2deg^. Removing one 
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Fig. 7. Covariance matrix correlation coefficients of w(0) from 
Jackknife estimates, for galaxy samples with 0.4 < z < 0.6 and 
Mg - 5 logh < -19.8. Left panel: red sample. Right panel: blue 
sample. 



sub-sample at a time for each Jackknife realisation, we compute 
the covariance matrix as 



M 1 '^ 

C(wi, Wj) = ^^^ Zj^"^'' ~ ^''^"^K' ~ ^J^ ' 



A^ 



(9) 



where A^ is the total number of subsamples, w the mean correla- 
tion function, and w' the estimate from the Z'*' Jackknife sample. 
Since w{9) is computed from a combined catalogue including all 
four fields, our Jackknife estimate leads to a noisy but fair esti- 
mate of cosm i c vari ance. We apply the correction factor given by 
iHartlap et al.l (l2007h to the inverse covariance matrix to compen- 
sate for the bias introduced by the noise. We show the correlation 
coefficient of the covariance matrix, r,j - Cjjl sJCaCjj for two 
red and blue samples, respectively, in Fig.|2l 

4. Filling haloes with galaxies: the halo occupation 
model 

4.1. The halo occupation distribution 

In order to relate galaxies to the dark matter haloes which host 
them we have implemented an analyti c model of galaxy clus - 
tering, the halo model (for a review see lCoorav & Shethll2002h . 
which contains at its core a prescription for how galaxies popu- 
late haloes, namely the halo occupation distribution. Our model 
follows closely the approaches used in recent works; further de- 
tails and references can be found in Appendix |A] The key as- 
sumption underlying our halo occupation distribution function 
is that the number of galaxies A^ in a given dark matter halo de- 
pends only on the halo mass M; it does not depend on environ- 
ment of formation history of the haloes. Furthermore, following 
IZheng et aP (l2005h . we express N(M) as the sum of two terms, 
corresponding to the contribution from the central galaxy N^ and 
the satellite galaxies N^'. 



N{M)^N^(M)x[\+N,(M)] , 
where 



Ndm 



Ns(M) 



1 



1 +erf 



log M - log M„ 

0"logM 



M, 



(10) 

(11) 
(12) 



The smooth transition for central galax ies expresses the u ncer- 
tainties in the galaxy formation process (IZheng et al.ll2007h . The 
factor Nc{M) for the satellite number in N(M) accounts for the 
fact that a halo cannot be populated by satellite galaxies without 
the presence of a central galaxy. 

For a given cosmology and dark matter halo profile, our 
model has five adjustable parameters. Mmi,, is the mass scale 
for which 50% of haloes host a galaxy. To reflect the scatter in 
the luminosity- halo mass re l ation, a smooth transition of width 
ciogM is used. IZheng et aP (l2007h show that if this scatter is 
small, the above expression takes the identical form as the distri- 
bution of central galaxies given a halo mass M, integrated over 
the entire luminosity range above the luminosity threshold. Thus 
Mmin also represents the halo mass scale for central galaxies, 
whose mean luminosity (Lc) is equal to the luminosity threshold 

This simple relation between Mmin and L^ is based on the 
hypothesis that stellar mass (or luminosity) has a power law de- 
pendence o n halo mass, which may not be exact over the entire 
mass range. iLeauthaud et al.l ( 1201 lal) recently proposed a model 
in which this relation assumes a more realistic form. The authors 
showed that in their model M^in and criog^ take different values 
than those computed with other models. However, in the mass 
range over which we compare our results with theirs in Sect. |6] 
(Mh ~ 10'-/2"'Mo), Mn,in values do not differ by more than 10% 
(see Fig. 3 in Leauthaud et al. 2011a). 

The number of satellite galaxies as function of halo mass 
follows a power law with slope a and amplitude Mi . Mi then 
represents the characteristic scale for haloes hosting one satellite 
galaxy. At lower masses, the dependence becomes steeper and 
the transition mass scale occurs at M ~ Mq. 

We show in Fig.[8]an example of measured w{0) and its best- 
fitting model, together with the best-fitting HOD function N{M). 
The total galaxy correlation function is the sum of two terms. At 
distances much smaller than the virial radius, the one-halo term 
contains contributions from galaxy pairs within a single halo, 
whereas at large distances the two-halo term contains contribu- 
tions from pairs in separate haloes. 



4.2. Deduced parameters 

From the HOD model we obtain deduced parameters describing 
galaxy properties. The mean galaxy bias bg at redshift z is the 
mass-integral over the halo bias bh (see Sect.|A]i weighted by the 
number of galaxies, 

r N{M) 



bJz) 



bt,(M,z)n(M,z)- 



(13) 



ngdz) 

The dark-matter mass function n is given in Eq. IA.6I and the halo 
bias bh in Eq. lA.14l The total number of galaxies is 



«gal(z) = I 



N(M)n(M,z)dM. 



(14) 



Similarly to the galaxy bias, the mean halo mass for a galaxy 
population is 

N{M) 



I 



<Mhaio)(z)= dMMniM,z) 



ngdz) 



The fraction of central galaxies per halo is 

NAM) 



I dMn 



(M,z)- 



ngai(z) 
Consequently, the fraction of satellite galaxies is 

fs(z) = l-Mz). 



(15) 



(16) 



(17) 
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Fig, 8. Example of a measured w{6) (for all galaxies in the redshift range 0.4 < z < 0.6 and for Mg — 5 log h < -19.8), as well as 
the best-fitting model, as described in sec.|4] Left: w{6) measurement and model. Right: N{M), showing the central term Ncem and 
the satellite term A^s- 



4.3. Population Monte Carlo sampling 

We use the "Populati on Monte Carlo" (PMC) technique to sam- 
ple li kehhood space JWraithet al.l 120091: iKilbingeret ai]|2010i 
1201 ih and have implemented it in a publicly-available code, 
CosMoPMCfl Contrary to the widely-used Monte Carlo Markov 
Chain method (MC MC), PMC is an adap tive importance - 
sampling technique (ICappe et al.l |2004 l2008h . which has two 
principal advantages: first, the parallel sampling algorithm com- 
bined with our fast HOD code means that each run can be 
quickly computed. Secondly, PMC does not have issues with 
chain convergence. Instead, the perplexity (defined below) is a 
diagnostic for the reliability of a given sampling run. 

Points are sampled from a simple, so-called importance sam- 
pling function q. Each sample point 6,, is attributed a weight w„ 
which is the ratio of the posterior n to the importance function. 



q(o„y 



J]^n 



1. 



(18) 



The initial importance sampling function ^ is a mixture of seven 
multi-variate Gaussians. A PMC run consists of a number of iter- 
ations, in the course of which the importance function is adapted 
to better match the posterior distribution. We run PMC for typi- 
cally 10 iterations using 10000 sample points per iteration. 

As a stopping criterion for the iterations we take the perplex- 
ity p which is defined as 



1 

N 



exp 



■ ^ w„ log 



(19) 



The perplexity p e [0; 1] is a measure of the distance be- 
tween the posterior and the importance function, and it ap- 
proaches unity if the two distributions are identical. The perplex- 
ity is a measure of the adequacy and efficiency of the sampling. 



Values of p above about 0.6 indicate a good agreement between 
the posterior and the importance function. 

The evaluation of the likelihood function, which is the 
most time-consuming process for most sampling tasks, is eas- 
ily parallelisable in importance sampling with little overhead. 
Furthermore, we have optimized several aspects of the numeri- 
cal computa tions of the HO D model. For example, we employ 
the FFTLog (lHamiltonl2000l) algorithm to perform Fourier trans- 
forms and make use of tabulated values to improve efficiency. 
The computation of the angular correlation function w on a range 
of scales is performed in under a second on a standard desk- 
top. This module is part of the latest public version (vl.l) of 
CosMoPMC. Typically, on an eight-core desktop (16 threads) a 
single PMC run takes around 90 minutes for a total of 100, 000 
sample points. Our PMC technique allows us to efficiently sam- 
ple the parameter space for the large number of galaxy samples 
used in this work. 



4.4. Likelihood function 

For each galaxy sample we simultaneously fit both the projected 
angular correlation function w and the number density of galax- 
ies «gai, by summing both contributions in log-likelihood: 

/ = Y,[^"'''m-w'^°'^\6d]{c-')..[w°'''{dj)-w"'°'''\0j)] 

["gal' - «"af '']^ 



'•] 



http : //cosmopmc . info 



(20) 



where n™^^"^ is given by Eq.[T4l at the mean redshift of the sam- 
ple. The data covariance matrix C is approximated by Eq.|9] The 
error on the galaxy number density cr„ _j contains Poisson noise 
and cosmic variance. The latter is estimated from the field-to- 
field variance between the four Wide patches. 
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The likelihood function is 
£ = exp(-;r'/2). 



(21) 



The product of prior P and likelihood defines the posterior 
distribution, n - P£.. We sample a five-dimensional space 
where {logi,)Mn,in,log,(,Mi,logioMo,a,o-iogM) is the parame- 
ter vector We use flat priors for all parameters; the ranges 
are logioMmin e [11; 15],logioMi e [12; 17],logioMo e 
[8; 15], a e [0.6; 2] and o-\o^m g [0.1; 0.6]. 

5. Results 

5. 1 . Clustering measurements 

Our measurements are presented in Figure |9] [10] and [TT] for 
the full, red and blue samples, respectively. For clarity, we 
show three samples selected with the same luminosity thresh- 
old {Mg - 5 logh < -19.8) and three samples covering the same 
redshift range (0.4 < z < 0.6). At the same luminosity and red- 
shift, red samples are more clustered on small scales than the 
full sample, and blue galaxies are the least clustered. At a con- 
stant luminosity threshold, w{6) increases with redshift for the 
full, red and blue samples. 

On the left panels of Figs. |9]and[T0] we plot our best fit w{d) 
models. Model parameters and deduced quantities are shown 
in Tables IB. II (full) and IB. 21 (red). An example of sample dis- 
tributions of best-fitting HOD parameters is given in Fig. IC.ll 
Since the halo model aims to reproduce the clustering in thresh- 
old samples as a sum of contributions from central and satellite 
galaxies, both types of galaxy must be present. We presume that 
in general, central galaxies are brighter than their satellites, so 
that the latter assumption is valid for full samples selected by lu- 
minosity threshold. Also, we assume that red satellites are sys- 
tematically associated with a (brighter) red central galaxy. On 
the contrary blue satellites may belong to a red central galaxy. 
For this reason, we defer fitting our model to our pure blue sam- 
ples. We plan to address the blue population in a future work, 
perhaps considering separat e halo occup ation functions for both 
red a nd blue populations CSi mon et al.l |2009; Skibba & ShethI 
I2009h . 

In most cases, we find a good agreement between the data 
and our model, with;^'^/dof from 1.77 to 10.21 for full samples 
and from 1.20 to 6.21 for red samples. We report a few large 
X^ values in some full samples, mostly at high redshift and for 
bright luminosity thresholds, where the galaxy number density 
is small. We suspect systematic photometric redshift errors to be 
responsible for an underestimated clustering signal and system- 
atic errors in the galaxy number density estimates (see Sec 13.2b . 
Perplexity (see Sect. 14.3b is larger than 0.6 for most samples, 
indicating that our model accurately describes the observations. 

5.2. Halo model fitting 

In the right panels of Figs.l9]and[T0lwe show our best fit N{M), 
the halo occupation function. First, we note that at a given red- 
shift interval, M^in, the average mass for which 50% of haloes 
contain one galaxy increases with luminosity. Secondly, we note 
that the evolution in N{M) with redshift at a fixed luminosity 
is less pronounced than the evolution with luminosity at a fixed 
redshift, whereas the amplitude of w{6) significantly decreases 
with redshift. 

Halo masses Mi and M^i„ are tightly constrained, whereas 
Mo, o-\ogM and a are poorly constrained. M^m and M[ mea- 
surements are displayed in the top panels of Fig. [12] with re- 



spect to luminosity threshold (Mg). For every sample, halo 
masses increase with luminosity. M^in covers the range 10" to 
10'** /i^'Mq. Ml is significantly larger than Mmm, ranging from 
]^q12.5 J.Q 1O'^/z"'M0. T hese trends are consisten t with those 
found in both th e local (TZehavi et aDl2010i l2005h an d distant 
Universe studies dZheng et al.ll2007l: I Abbas et al.ll20Toh . At high 
luminosity thresholds, where the number of objects is small, 
error bars on our best-fitting parameters are correspondingly 
larger. This trend is especially pronounced for a, where the small 
number of satellites in high mass haloes makes it diflicult to 
constrain this parameter However, we note a increases with lu- 
minosity threshold, suggesting that for massive haloes, bright 
galaxies will be more numerous than faint ones. criogM remains 
relatively stable for red galaxies (~ 0.3), as also seen in other 
data sets and numerical simulations. For the full sample, val- 
ues are slightly larger (up to ~ 0.5 in some cases). Our fitted 
o"iogM values are slightly larger than previous works. We also 
made fits where values o-\ogM > 0.6 were allowed, and in some 
cases large values (> 1.0) were found. We therefore decided to 
impose the restriction that a-\ogM < 0.6. Sample incompleteness 
(due to errors in the photometric redshifts), leading to missing 
central galaxies, could explain such high o-io„m values. 

At constant luminosity threshold, in both red and full sam- 
ples. Ml and M^in decrease with increasing redshift. This effect 
is partially due to a selection effect caused by the dimming and 
reddening of stellar populations with time; for a constant lumi- 
nosity cut one selects less and less massive galaxies at higher 
redshifts. We will return to this measurement in Sect. [6] where 
we attempt to separate this passive stellar evolution effect from 
the intrinsic stellar-to-halo mass evolution. 

At faint luminosities, Mmi,, values are higher for red galaxies 
than for full samples. As we explained previously, one assump- 
tion of our model is that our red threshold samples have lumi- 
nous red central galaxies. On the other hand, it is also possible, 
especially for fainter samples, that blue satellite galaxies inhabit 
haloes with red central galaxies. In such a case for the faint red 
galaxy sample, a red central galaxy will be seen "alone" in a 
more massive halo. This would also increase the halo mass esti- 
mate found for red central galaxies as compared to those in the 
full sample. 

The bottom panels of Fig. [12] show Mi with respect to M^i^ 
for full and red samples. The parameters are strongly correlated 
and we fitted a power law (dashed line) in each redshift bin. For 
the full sample, we compare our results with a linear relation. 
Ml = 17 X Mmin (dotted line), based on the c oefficient observed 
in the SDSS [z ~ 0.1) bv lZehavi et"an(l2010l) . Our results agree 
with them in the lowest redshift bin, z ~ 0.3, but depart from this 
value at higher redshift. 

We find that although Mi /Mmin relation does not vary signif- 
icantly with luminosity threshold it does show some indications 
of evolution with redshift. We note that Mi /Mmin decreases with 
redshift for M„iin < IO'^'/i^'Mq and increases with redshift at 
higher mass. For red galaxies. Mi is no longer proportional to 
Mmin but follows a power law. For M,nin ~ IO'^/j^'Mq, the ratio 
Ml /Mmin is ~ 12 (as compar ed to ~ 10 for luminous red galaxies 
found by IZheng et aDl2009h and shows a much smaller redshift 
evolution over the entire mass range. 

Figure [13] shows the relationship between Mmin and the 
galaxy number density «gai for red and full samples. As one 
would expect, progressively rarer objects are found in more mas- 
sive haloes. We do not detect a significant evolution of the slope 
of our Mmin versus n„ai relationship wit h redshift. For red g alax- 
ies, we overplotthe measurements from lBrown et al.l ( |2008|) with 
which we are in excellent agreement, even though it covers 
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Fig. 9. w(&) measurements (left) with their best fit models for all galaxies, and the resulting galaxy occupation function (right) with 
respect to halo mass. We illustrate the redshift evolution (top) in the luminosity threshold M^ - 5 log/i < -19.8, and the luminosity 
dependence (bottom) in the redshift range 0.4 < z < 0.6. 



a much smaller area. For t he full galaxy sarn ple we compare 
with SDSS measurements (IZehavi et al.ll2010l) . The agreement 
is good at high surface densities; however we note a departure 
for low density samples. This difference suggests that at con- 
stant density, galaxies tend to reside in more massive haloes at 
low redshift (z ~ 0) than at higher redshift (z > 0.6). Part of 
this diff'erence could be caused by systematic photometric red- 
shift errors. As we discussed in Sec 13.21 bright samples may 
be contaminated by fainter objects, and the combined effect of 
higher density and lower clustering signal will result in reducing 
Mtnin. In the most pessimistic case outlined in Sec l3.2l where the 
density estimate would be overestimated by a factor 3, log Mmin 
would be biased low by a few sigmas (up to ~ 0.5). 



The relationship between Mmin and «gai can be well approx- 
imated by a power law and we fitted galaxy samples at all red- 
shifts simultaneously. We found: 



M± 



10'°° X \nfJ(}^^ Mpc-^^)]-°«4 /,-' Mq , 



for all galaxies and 



xlO.3 , 



red 
'eal 



3M-0.74 , -1 



M'^^ = 10'"-" X [n'^^JilfM^c-')]-"'"!!-' M( 



O^ 



for red galaxies. 



(22) 



(23) 



6. Discussion 

The key issue we would like to address is to understand how 
galaxy formation and evolution depend on the properties of the 
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Fig. 10. w(ff) measurements (left) with their best fit models for red galaxies, and the resulting galaxy occupation function (right) with 
respect to halo mass. We illustrate the redshift evolution (top) in the luminosity threshold Mg - 5 log/i < -19.8, and the luminosity 
dependence (bottom) in the redshift range 0.4 < z < 0.6. 



underlying dark matter haloes. Ideally, to relate halo masses 
to their stellar content, one would select our samples by stel- 
lar mass. However, stellar mass estimates computed from our 
five-band optical data alone would suffer from large uncertain- 
ties. Keeping this limitation in mind, we adopt an intermediate 
approach where we determine an empirical relation to convert 
the observed luminosity threshold of each sample to an approx- 
imately mass-selected sample using a reference sample with ac- 
curate stellar masses. This method allows to transform luminosi- 
ties into masses and provides an estimate of the likely errors in 
such a transformation. 



6.1. Transforming to stellar mass threshold samples 

As a consequence of passive stellar evolution, the characteristic 
rest-frame M^-band luminosity of a given galaxy t ype evolves 
signi ficantly over the redshift range 0.2 < z < 1.2 (lllbert et al.l 
i200llZucca et al.ll2006l:lFaber et al.ll2007h . 

In order to account for these effects, we have established 
a simple approximation to relate stellar m ass and luminosit y 
based on COSMOS 30-band photometry (lllbert et alJ l2OI0l) . 
We used the COSMOS B-band luminosities in order to be as 
close as possible to our M^-band selected samples. Red galax- 
ies are selected as M(NUV) - Mr > 3.5 and blue galaxies as 
M(NUV) - Mr < 3.5. We then fit the stellar mass-to-luminosity 
ratio with power-law functions as function of redshift between 
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z = and z = 1.5 for five stellar mass bins from 10^" to 10"'^. 
These results are displayed in Fig. [14] 

Due to the ageing of stellar populations, the average stel- 
lar mass-to-light ratio increases with time. For red galaxies, the 
slope depends weakly on stellar mass at intermediate masses 
(10''-^-10"-^/i"'Mo) and ranges from -0.56 to -0.44. The slope 
is steepest in the lowest mass bin, but the number of objects in 
this bin is small. Based on these results we adopt the following 
relationship for the red galaxy population: 



log 



M„ 



■^B /red 



= log 



M„ 



-(0.5±0.1)z, 



(24) 



B /i-ed„'=0 



where the 0.1 error accounts for the scatter between mass bins. 
This evolution is consistent with the luminosity function evolu- 
tion me asured for red galaxies in the DEEP2 and COMBO 17 
surveys (iFaber et al .1120071: IWillmer et al.ll2006h . where the char- 
acteristic absolute magnitude M^ decreases by about 1.2-1.3 
mag per unit redshift (Mg cc -2.5 log Lg). 

For the blue sample, the evolution in redshift is more 
pronoun ced than for the red sample, as also seen i n 
DEEP2 (Wi llmer et all |2006|) and VVDS dZucca et alj |2006|) . 
Unfortunately, in this case, the slope depends on stellar mass, 
ranging from -0.69 to -1.10 in the mass range [10^'^, lO"'^], 
and has a larger scatter The full sample is dominated by blue 
galaxies at faint luminosities but by red galaxies at bright lumi- 
nosity (see Fig. [3]). Therefore a simple relation between luminos- 
ity and stellar mass valid for all luminosities cannot be obtained. 
For the full sample we choose to apply the same correction as for 
the red sample. As blue galaxies evolve more rapidly than red 
ones such a correction underestimates the expected evolution of 
the full samples at faint luminosity and therefore cannot be used 
to construct stellar-mass selected samples. In the following sec- 
tions, we discuss how our results depend on this correction. 

In the following sections we use these expressions to cal- 
culate "corrected" luminosities, denoted as U. Note that these 
corrections are equivalent to shifting observed galaxy luminosi- 
ties to redshift z - 0. Our results are displayed as a function 



of relative luminosity threshold L' /L^, where L, corresponds to 
M* - 51og/! = -19.81 and -19.75 for the full a nd red sam- 
ples, respectively, measured in the local Universe dFaber et al.l 
l2007l and references therein). It is also worth mentioning that if 
galaxies would experience a pure passive luminosity evolution 
with cosmic time, this would result in constant HOD parameters 
for a given L' /L^ as function of redshift. However, deduced pa- 
rameters (galaxy bias, mean halo mass and satellite fraction) are 
expected to vary due to the evolving halo mass function and halo 
bias (see Eqs. [T3l - [T7l i. 



6.2. The stellar-to-halo mass relationship 

In our model, M^in represents the hosting halo mass for central 
galaxies whose median luminosity equals the luminosity thresh- 
old of the sample (see Sect. 14.1b . The ratio L' /M^^ is therefore 
the ratio of luminosity to dark matter mass for central galax- 
ies. IZehavi et al.l (1201 Ol) showed that the relation between cen- 
tral galaxy luminosity L^ (corresponding to our L') and the halo 
mass could be approximated by: 






(25) 



where A, om and Mi are free parameters and Mh represents M„im 
01 Ml. 

This expression encapsulates the idea that there exists a 
"transition halo mass" (~ M,) where the ratio of central galaxy 
luminosity to dark matter halo mass reaches a maximum. This 
transition mass represents the halo mass for which baryons have 
been most efficiently converted into stars until the time of ob- 
servation, resulting from past and on-going star formation over 
cosmic history. 

We fit Eq.|25]for full and red samples for each redshift slice. 
For our two highest redshift bins where the number of points are 
less than the number of free parameters, we perform a simulta- 
neous fit over both redshift bins. Results are plotted in Fig. [15] 
and best-fitting parameters given in Table [j 
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In our lowest redshift bin (0.2 < z < 0.4), our full sample 
has best-fitting parameters slightly di fferent than in the SDSS 
dZheng et al.ll2007i; IZehavi et al.ll2010l) . A = 0.20 + 0.01 (com- 
pared to 0.32) and M^ = 2.24 + 0.1 x 10"/r'Mo (compared 
to 3.08 X 10"). We note that the difference is likely due to the 
different selection {Mr in the SDSS). At high luminosities, the 
difference with their relation increases. For a given luminos- 
ity, we find lower halo masses, L^ oc ^o.5±o.oi^ compared to 
Lc oc ^0 26-0.28^ which follows the trends observed in Fig.[T3l 
and might be due to contaminated samples at low d ensity. In the 
red sa mple, we measure L^ oc ^o.3i±o.i6^ whereas [Zheng et all 
(120091) found L^ oc m"^ for luminous red galaxies in the SDSS. 
Both Mmin and M\ increase with redshift for the full sample. The 
dependence with redshift is weak for red galaxies, although we 
detect a slight increase in the very highest redshift bins. 

For the full sample we assumed a correction based on the 
evolution of the stellar mass-to-luminosity ratio of red galax- 
ies (Eq. l24b . As the characteristic luminosity of blue galaxies 
evolves more rapidly than red ones and as the faint sample is 



dominated by blue galaxies a legitimate concern is the depen- 
dence of our results on our corrections to stellar mass. Several 
corrections assuming the blue galaxy luminosity evolution (us- 
ing slopes given in Fig. [T4l) were tested. If a larger correction 
is applied, M\ and Mmi„ show a similar but more pronounced 
trend. For this reason we conclude that the increase of halo mass 
parameters with redshift for the full sample is a robust result (al- 
though the exact amount of evolution observed depends on the 
correction we apply). 

6.3. Redshift evolution of the SHMR and comparison with 
iiterature vaiues 

To study the redshift evolution of the stellar-to-halo mass ratio, 
we consider the relationship between corrected galaxy luminos- 
ity and the host halo mass from Eq. |25) 



Mh M,\MJ ^\ Mh 



+ 1 



(26) 
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where the maximum is located at: 



Ml 



peak 



M, 



1 -Om 



(27) 



The right panels of Fig. [15] show the L'^jM^ia relation 
for full and red samples as function of halo mass fitted 
by Eq. |25] Best fitting parameters are listed in Table IC.6I 
For less massive haloes, the relative stellar mass content 
increases with halo mass and reaches a maximum around 



M, 



\Q h Mq. In more massive haloes, the stellar- 



to-halo mass ratio sharply decreases. Our measurements are 
consistent with recent observations measu ring the stellar-to- 
halo mass ratio as function of halo mass (iMandelbaum et all 



l2006l:IZheng et alJ2007HFoucaud et alJ2010t]zehavi et alJ2010l: 
More et al..,201 iLIWake et al.l201 itJLeauthaud et al.ll201 Ibl) . and 



"abundance-matching" techniques based on A^-body simulations 
(IConrov & WechsleJ l2009t JMoster et aP 120101: iBehroozi et al.l 



I2OIOI) . Simila r trends are found in semi-analytic simulations 
(see, for e.g., ISomerville et alj r2008h , in which physical pro- 
cesses (such as supernovae and AGN feedback) preventing gas 
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from cooling and quenching star formation must be taken into 
account to match observations. 

Our measured peak values M^'^^ (Eq. |27] with error bars 
from the nearest data points) are plotted as a function of red- 
shift in Fig. [16] Comparing to literature measurements at lower 
redshifts, the peak stellar-to-halo mass relation is M^'^ - 
4.54 X 10"/i ~'Mn, similar to th e local Universe value (^ 4.2 x 
10"/!"'Mh, IZehavi et alll2010l) . For our full sample, the peak 
location gradually increases with redshift. 

At h igher redshif t s, we com pare our meas u remen ts with 
those of lMosteret"an (1201 Oh and (Behroozi et al.l (l2010h . These 
authors derived the stellar-to-halo mass relation using an 
abundance-matching technique, linking dark matter halo mass 
functions from A^-body simulations with observed galaxy stellar 
mass functions. JMoster et al.l (1201 Ol) estimate the redshift evolu- 
tion of the stellar-to-halo mass ratio (their Table 7) and from this 
we derived analytically the position of the peak, as well as asso- 
ciated errors. Values fromJBehroozi et al. (2010) at redshift 0.1, 
0.5 and 1 .0 were supplied to us (P. Behroozi, private communica- 
tion). Both lMoster et al.l and lBehroozi et aLT s measurements are 
based primarily on galaxy stellar mass measurements covering 
relatively small areas (a few hundred square arcminutes). Our 
full sample is in excell ent agreement withiBeh roozi et al. but dif- 
fers significantly from JMoster et al.L although in both cases the 
two data sets show the same general trend: Mj^*^" shifts gradually 
towards higher halo masses at higher redshift. The Moster et al. 
points do not seem to agree very well with most observations and 
moreover overestimate the position of the peak at low redshift. 

Measurements in t he COSMOS field are also shown 
(JLeauthaud et al.ll201 Ibl) : their points are in general larger than 
most measurements, with the exception of Moster et al. . It has 
already been not ed that there are several rich structures in the 
COSMOS field (iMcCracken et al.l l2007t iMeneux et al.l 120091: 



Ide la Torre et al.l 1201 Oh which lead to higher correlation func- 
tions, especially at z ~ 0.8; cosmic variance could perhaps be 
the origin of this discrepancy. 

For red galaxies, the peak is at larger halo masses, Mj^*^^ = 
20.6 X 10"/2"' Mq. This is not surprising as Mmi,, is larger for red 
galaxies at faint luminosities. The red galaxy peak position un- 
dergoes only little evolution to z ~ 1 . Our intermediate-redshift 
red ga laxy points are co nsistent with higher redshift measure- 
ments ( IWake et al.ll2011l) : these authors used a near-infrared sur- 
vey to select galaxies by stellar mass, and we would expect their 
selection to be dominated by red galaxies. 

The results presented here represent the first time a single 
data set covering a statistically significant area has been used to 
derive the evolution of the peak position as a function of red- 
shift and galaxy type in a self-consistent manner. The peak po- 
sition, MJJ*^^ , can be interpreted as representing the halo mass at 
which the stellar mass content is most efficiently accumulated in 
haloes, either by star formation or merger processes. The move- 
ment of this peak towards higher halo masses at higher redshifts 
is consistent with a picture in which the sites of efficient stel- 
lar mass growth migrate from low-mass to high-mass haloes at 
higher redshift. This "anti-hierarchical" evolution, frequently re- 
ferred to as "ha lo downsizing", res embles closely the scenario 
first sketched bv lCowie et al.l ( Il996l) . in which the maximum lu- 
minosity of galaxies undergoing star-formation declines steadily 
to the current day. Conversely, the fact that we observe constant 
transition masses for the red population suggests that the stellar- 
to-halo mass ratio does not significantly evolve in haloes with 
mass ~ IO'^/i^'Mq and supports observations that intermediate- 
mass red galaxies have experienced no significant stellar mass 
growth by star formation since z ~ 1.2. 

We tried several corrections to convert to stellar mass thresh- 
olds based on the slope measured for blue galaxy evolution 
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(see Fig. [T4l) . and we found that the peak luminosity depends 
only weakly on these corrections. Only the highest-redshift point 
changes by around one sigma if we apply the extreme "blue 
galaxy" correction to our full sample. For this reason we con- 
clude that the increase of MJJ*^" with redshift for the full sample 
is a robust result. 



6.4. Redshift evolution oflialo satellite fraction 

The fraction of galaxies which are satellites in a given dark 
matter halo is a sensitive probe of the past evolutionary his- 
tory of the halo which may be modified by processes such 
as major mergers (which can decrease the satellite frac- 
tion with time) or environmental effects which operate pri- 
marily to reduce the number of satellites, such as, but not 
limited to, the res triction of gas supply or " strangulation" 
( Larson et al.l Il98d) or ram-pressure stripping JGunn & GottI 
Il972h . Satellite galaxies may also merge with their central 
galaxy. Observational evidence based on pair fraction measure- 
ments and theoretical studies using numerical simulations now 
shows that mergers seem to play a significant (although not dom- 
inant) role in the evolution of massive galaxies since z ~ 1 



(White et al. 2007 ; Wang & Jing||2()10l:lvan Dokkum et al.ll2010l: 
fkobaina et al...2010: .Zehavi et al.ll201 11) . Our objective in this 
section is to see if measurements of the satellite fraction can pro- 
vide new insights into the evolutionary history of haloes. 

The left panels of Figure [TT] show the satellite fraction /,, 
computed using Eq. [17] As before, each sample is converted to 
an approximately mass-limited sample by "correcting" for pas- 
sive luminosity evolution (see Sect. 16.1b . At all redshifts, as the 
luminosity threshold increases, corresponding to higher mass 
thresholds, the satellite fractions become progressively smaller. 
The similarity between the full and red samples is simply a con- 
sequence of the dominance of red galaxies at bright luminosity 
thresholds. 

At fainter luminosity thresholds, the satellite fraction rises 
and eventually reaches a plateau. This effect can be understood 
by considering both the halo mass function and the halo occu- 
pation function: the increase of the latter with decreasing lumi- 
nosity threshold (see Fig.[8]i is rapidly compensated by a sharp 
decrease in massive halo number density, so that faint galaxies 
are preferentially central galaxies in smaller haloes rather than 
being satellites in more massive haloes. 

We compare our results with satellite fractions found in 
SDSS (z ~ 0.1) and DEEP2 (z ~ 1). To explore a complete 
range in redshift and luminosity, we extrapolate our best-fitting 
HOD parameters measured in the redshift range 0.4 < z < 0.6; 
that is, at a given luminosity, we keep the HOD parameters fixed 
and compute Eq. [T7]at different redshifts. The choice for this in- 
termediate redshift range is motivated by the fact that it contains 
the most robust photometric redshift estimates and our measured 
points span a large range in luminosity. Each HOD parameter set 
is constructed from Mmi,, and Mi interpolated in luminosity us- 
ing Eq. [25] with best-fitting parameters given in Table IC.6I For 
red galaxies, Mq is set to M^in/lOO, ctiq^m is fixed to 0.3 and 
a = 1 . 1 for L < L» or ff = 1 .4 otherwise. For all galaxies, Mq 
is set to Mmin/100, (TiooM = 0.4 and a = 1.1 for L < L^, and 
o"iogM - 0.3 and o- = 1.3 for brighter luminosity thresholds. 

In the middle panels of Fig.[T7]we compare our extrapolated 
satellite fractions for all galaxies (in dashed lines) with results 
derived from SDSS and D EEP2, using a HOD analysis similar 
to ours dZheng et alj|2007l) . Our satellite fractions are consistent 
with these measurements at both redshifts, although we slightly 



over-predict satellite fractions at z ~ 0.1 for samples brighter 
than ~ 0.5L,. An overall increase of a factor of two is measured 
between z ~ I and z ~ for samples with L > Lt,we note that 
this value is lower than the factor of three from z ~ 0.5 to z ~ 
measured in the VVDS d Abbas et al.ll2010h . 

We fit the HOD parameters at each redshift independently. 
Therefore, our model provides us with "snapshots" of the occu- 
pation function at various epochs. If any physical process had 
played a significant role in reducing the satellite number density 
between z ~ 1 and z ~ 0, our fixed-HOD model calibrated at 
z ~ 0.5 would over-predict the satellite fraction observed in the 
local Universe and under-predict the satellite fraction at redshift 
one. The excellent match found for bright luminosity thresh- 
old samples shows that our satellite fraction measurements fol- 
low closely a redshift-independent HOD. Since a redshift evo- 
lution of the merging rate would change the occupation func- 
tion, we therefore conclude that our measurements imply a con- 
stant merging rate since z ~ 1, which confirms the constant (al- 
beit minor) role of merging for inter mediate-mas s galax ies (e.g., 
lBundvetal.1120091) . As outHned by IWhite et all (|2007|) . if halo 
merging was not followed by galaxy merging, we would observe 
an increase of halo masses at constant galaxy density. Instead, no 
redshift evolution of the Mn,in-«gai relationship (see Fig. [13]) sup- 
ports a scenario with constant galaxy merging with time. 

We compare our red satellite fraction with the galaxy-galaxy 
lensing analysis from SDSS (Mandelbaum et al. 2006). These 
samples were constructed from luminosity slices, rather than the 
threshold samples we use here. We then computed the number of 
satellites per luminosity bin as the difference of satellite fractions 
between two threshold samples: 



fs 
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(28) 



As we found for the full galaxy sample, the good agreement sup- 
ports a non-evolving HOD from z = 0.5 to z = 0. 

Finally, in the right panels of Fig. [17] we provide extrap- 
olated satellite fractions as function of redshift for several lu- 
minosity thresholds. It is interesting to examine our results in 
the context of recent results studying t he role of satellite galax - 
ies in galaxy formation and evolution dPeng et al.ll2010i 1201 lb . 
These works show measurements from numerical simulations 
and SDSS group catalogues which indicate that the fraction of 
satellites at a fixed overdensity is independent of halo mass and 
redshift. In this context, the growth of the satellite fraction we 
observe is simply a consequence of average growth in overden- 
sity in rich structures between z ~ 1 to z ~ 0, reflecting the 
growing importance of environment at low redshift. 

6.5. Galaxy bias 

We compute the galaxy bias (Eq.[T3]) from our best-fitting HOD 
parameters. We fit the bias measurements as fun ction of lumi- 
nosity and redshift, adopting a linear relationship dNorberg et al.l 
,200lh: 



L' 

bg{> L) =abias + ^bias— , 



(29) 



where L' represents the corrected luminosity (computed from 
Eq.[24]for the full and red samples). We fit the two highest red- 
shift bins together (0.8 < z < 1.0, 1.0 < z < 1.2), as only one 
point is available in the highest one. Results are plotted in Fig.fTSJ 
and best-fitting parameters are given in Table lC.7l 
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Fig. 17. Satellite fractions as function of luminosity and redshift, for all and red samples (top and bottom panels). The left panels 
show the measured satellite fraction as function of corrected luminosity. Dashed lines in the middle panels show satellite fractions 
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ratios. 



This parameterised bias relation reproduces our results best 
at z < 0.6. At a fixed redshift, for L'/i* < 1, the bias is a 
very weak function of luminosity; for L'/i* > 1, bias depends 
strongly on luminosity. For the faint, low-redshift full-sample 
galaxies, dominated by more weakly-clustered bluer galaxies, 
b ~ I; the faint red samples have higher biases, of around 
b ~ 1 .5. These results are a consequence of the higher clustering 
amplitudes observed for red samples compared to the full galaxy 
sample for L'/L, < 1, and are consistent with similar measure- 
ments m ade by sm aller surveys at intermediate redshifts (see for 
example [Meneux e t al. 2009). 

The right panels of Fig.[T8]illustrate the redshift evolution of 
the bias at several luminosity thresholds. We compare our results 
(solid lines) derived from Eq. |29] with a galaxy bias evolution 
computed using a constant set of HOD parameters evaluated at 
z ~ 0.5 (see Sec l6.4b . The slight difference observed between the 
measured HOD and the HOD fixed model at z = 0.5 comes from 
the fact that we set cr\o^M to 0.3 and a to 1.0. For the full galaxy 
population, a model with fixed HOD parameters provides an ex- 
cellent fit to the observations, with more luminous objects un- 
dergoing a much more rapid evolution in bias. For the red galaxy 
population, there is some evidence that the bias evolution is more 
rapid than the model with fixed HOD parameters. 



6.6. Mean halo mass 

From Eq.[T5]we compute the average dark matter halo mass for 
each of our samples. Note that this quantity represents the halo 
mass averaged over samples containing galaxies above the lu- 
minosity threshold and differs from M,nin which represents the 
halo mass for central galaxies whose luminosity corresponds to 
the luminosity threshold. As before, after correcting luminosities 
to a corresponding stellar mass threshold, we fit the mean halo 
mass as function of luminosity and redshift, adopting a linear 
relation: 



log<Mhalo>(> L') = Cfhalo + ^halo 



L' 



(30) 



where L' is the corrected luminosity. Once again, we fit the two 
highest redshift bins simultaneously. Results are displayed in 
Fig. [T9] and best-fitting parameters are given in Table lC.Tl 

For fainter samples with L' < 2L, the mean halo mass 
gradually increases with luminosity and changes more rapidly 
for brighter samples. Faint red galaxy samples on average have 
higher halo masses than the full sample which is dominated by 
bluer galaxies at these luminosities. This is consistent with a sce- 
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nario in which high numbers of red satellite galaxies in faint 
samples reside in more massive haloes. 

In the right panels of Fig. [19] we show the redshift evolution 
of the mean halo mass for several luminosity thresholds. As be- 
fore, we compare mean halo masses derived from our measured 
HOD parameters (solid lines) with a model of fixed halo model 
parameters set measured at redshift 0.5. At all luminosities, the 
mean halo mass decreases with increasing redshift. This is an in- 
dication of the fact that haloes merge through cosmic history and 
the effective mass of haloes grows with time. In bright samples, 
the mean halo mass evolves more slowly than in fainter ones. At 
low luminosity thresholds and high redshifts, both the full sam- 
ple and red galaxy sample deviate slightly from our model with 
fixed HOD parameters. 



7. Summary and conclusions 

We have made one of the most precise measurements to date of 
the angular correlation function w and its dependence on lumi- 
nosity and rest-frame colour from z ~ to z ~ 1. Our mea- 
surements were constructed from a series of /' < 22.5 volume- 



limited samples, containing ~ 3 x 10* galaxies in four indepen- 
dent fields. These samples were created using accurate five-band 
photometric redshifts in the CFHTLS Wide survey calibrated 
with ~ 10"* spectroscopic redshifts. A cross-correlation analysis 
showed consistency with a relatively small contamination frac- 
tion arising from scatter between photometric redshift bins. 



We interpreted these measurements in the framework of a 
model in which the number of galaxies inside each dark mat- 
ter halo is parametrised by a simple analytic function (the halo 
occupation distribution) with five adjustable parameters. In this 
model, central and satellite galaxies are treated separately. For 
each galaxy sample, we performed a full likelihood analysis and 
explored model parameter space using an efficient Population 
Monte Carlo analysis. The very large survey volume probed by 
the CFHTLS Wide allowed us to place robust constraints on 
the redshift evolution of the halo model parameters, and make 
accurate estimates for halo properties of many different galaxy 
samples over more than two orders of magnitude in luminos- 
ity and three orders of magnitude in halo mass, from 10'^ to 
lOi5(/riM0). 
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Using 30-band photometric data from the COSMOS survey, 
we derived an empirical relation between stellar mass, luminos- 
ity and redshift and used this to convert our luminosity-selected 
samples to stellar-mass limited samples. All of our conclusions 
(following below) are independent of this correction. 

These are our principal results: 

1 . For a given luminosity threshold, galaxies with redder rest- 
frame colours are more clustered than bluer ones. Clustering 
strength increases with increasing luminosity, reflecting that 
bright galaxies reside in more strongly clustered massive 
haloes. 

2. We consider the redshift and luminosity evolution of halo 
parameters Mi , representing the characteristic halo mass re- 
quired to host at least one satellite galaxy, and Mn,in, repre- 
senting the halo mass where on average 50% of haloes con- 
tain one galaxy. For our full sample. Mi is closely approxi- 
mated by ~ 17 X M,nin. This ratio becomes smaller at higher 
redshifts for lower-mass haloes. For red galaxies, Mi/M^i„ 
is ~ 12 at intermediate mass scale, and remains constant for 
all halo masses over 0.2 < z < 0.8. This lower ratio for 



redder galaxies reflects the higher abundance of satellites in 
these haloes; in general, redder populations reside in more 
massive haloes. 

By fitting a simple analytic relation between central galaxy 
luminosity and halo mass to our observations, we find a max- 
imum in the stellar-to-halo mass ratio. In the lowest redshift 
bin z ~ 0.3, this peak mass is M^'^ = 4.5 x IQ^^H^^Mq 

for the full samples, and M^'"'^ = 21 x IGI'/j-'Mq for the 
red samples, respectively. This transition represents the halo 
masses where baryons were most efficiently converted into 
stars, and is in good agreement with measurements from 
other studies. 

For the full sample, M^'^^ shifts to higher halo masses at 
higher redshifts. For the red galaxy sample, the peak po- 
sition evolves less rapidly with redshift. These results can 
be understood qualitatively from the lack of on-going star- 
formation in the red galaxy population which means that the 
stellar mass content in these massive haloes changes very 
slowly. For the full galaxy sample, which is expected to con- 
tain galaxies still undergoing active star-formation, the shift 
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of the transition mass to progressively more massive haloes 
at higher redshifts is a manifestation of "anti-hierarchical" 
galaxy formation. These results also indicate that the mas- 
sive, passive galaxies in our survey are already fully formed 
byz~ 1.2. 

5. At increasing luminosities, the fraction of satellite galaxies 
rapidly drops to zero. For less luminous galaxies, the satellite 
fraction is higher in our red sample than in our full sample, 
reflecting the high number of faint red s atellites, consistent 
with observations in the local Universe (iMandelbaum et all 
l2006l:IZehavi et alll2010h . For our full sample, the number of 
satellites increases by a factor of two from z ~ 1 to z ~ 0, 
consist ent with a combine d study of galaxies in DEEP2 and 
SDSS dZheng et alj|2007h . This result suggests that the oc- 
cupation function for satellite galaxies inside intermediate- 
mass haloes (~ lO'^/i^'M©) remains constant with cosmic 
time, and implies that the galaxy merging rate does not 
change since z = 1 . 

6. Above the characteristic sample luminosity, the mean galaxy 
bias bg increases rapidly. Faint red samples have bg ~ 1.5, 
compared to Zjg ~ 1 for the full sample. Redder samples also 
have higher average halo masses, which decrease towards 
higher redshifts. 

7. We compare the measured evolution of our mean halo 
masses, biases and satellite fractions with a model assum- 
ing constant HOD parameters over our redshift range. This 
model provides a remarkably good fit to our measurements 
for massive haloes. We note that lower mean halo masses 
(particularly in the full sample) show a slower decrease com- 
pared to a redshift-independent HOD model. 



Systematic photometric redshift errors have little impact on 
our faint luminosity-threshold samples, for which we find no ev- 
idence for incompleteness larger than the field-to-field variance 
estimates. Our very brightest samples however may be affected 
by incompleteness or contamination from fainter objects. We 
have shown that these errors could reduce our measured clus- 
tering amplitudes and consequently could lead to an underes- 
timation of halo masses. The effect will be more important in 
contaminated samples, where the higher galaxy number density 
will lead to lower halo masses. In the worst case, log Mmin might 
be underestimated by ~ 0.5. 

Because of the lack of deep near-infrared data in the 
CFHTLS, we cannot calculate reliable stellar masses and cleanly 
separate the active and passive populations in mass-selected 
samples. Such a data set, perhaps combined with approximate 
local density estimators calibrated using spectroscopic redshifts, 
and a refined halo model would allow us to understand in greater 
detail the physical origin of the shift of the transition masses to 
higher halo masses observed here. Furthermore, given the fact 
that the peak of the star-formation and mass assembly takes 
place in the redshift range 1 < z < 2, it is important to push these 
studies to higher redshifts by constructing larger, wide-area sur- 
veys with deep near-infrared data. 

In our model, we have assumed a halo mass function, halo 
density profile, and a cosmological model, although it is well 
known that the evolution of galaxy clustering depends on these 
parameters. Combining galaxy-galaxy lensing, galaxy clustering 
and mass-function measurements from an expanded CFHTLS 
including near-infrared data should provide new insights on both 
cosmology and galaxy evolution. 
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Appendix A: The halo model 

The halo occupation distribution (HOD) model 
("Berlind & Weinberg! l200l iKravtsov et all l2004t IZheng et al.l 
12005) allows one to derive physical information about galaxy 
populations and the dark matter haloes which host them. The 
HOD prescription is based on the halo model, which describes 
how dark matter is distributed in space. In this framework, all of 
the matter is assumed to reside in virialised haloes. The HOD 
parametrisation specifies how many galaxies populate haloes, 
on average, as function of halo mass. Accordingly, the number 
of galaxies per halo A^ only depends on the mass M of the halo. 



A.1. Dark-matter halo model 

The three ingredients to the halo model of dark matter are 
the halo mass function, the halo profile and the halo bias. 
The dark matter halo a bundance can be inferred using the 
iPress & Schechten (Il974h approach where dark matter collapses 
into overdense regions above the critical density Sc, linearly 



evolved to z = 0. The mass function, which is the halo number 
density per unit mass, can be parametrized as 



«(M,z)dM= — /(v)dv, 



(A.l) 



where po is the mean density of matter at the present day. The 
new mass variable v writes 



Sciz) 



D(z)o-{M) 



(A.2) 



and characterizes the peak heights of the density field as function 
of mass andredshift. The linear critical density 6c depends on the 
adopted cosmology and redshift; we use the fit ting formula from 
[Weinberg & K amionkowskil (l2003h : see also iKitavama & Sutol 
dl996D, ' 
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(127r)^/^[l + 0.013 logioQm(z)]; 



f^m(z) = nm(l+z)' 



Ho 



H(z) 



(A.3) 
(A.4) 



D{z) is the linear growth factor at redshift z, and cr(M) is the 
rms of density fluctuations in a top-hat filter of width R = 
{3M/4-7rpo)^^^, computed from linear theory. 



cr^iM) 






^''"^'''^wHkR), 



27r2 



(A.5) 



where W(x) = (3/jc^)[sin jc- x cos x]. For the mass function f(v), 
we choose the parameterisationbv lSheth & Tormenld 19991) . cal- 
ibrated on simulations: 



yf(y) 



:Aj^[l+(flyV]exp 



av 



(A.6) 



(A.7) 



where the normalisation A is fixed by imposing: 
jn{M,z)=dM = j f(v)dv = 1 , 

We adopt the values p - 0.3 and a - 1/ y2. If not indicated 
otherwise, all integrals over the mass function are performed be- 
tween Mio„ = lOVi-'Mo andMhigh = lO'^/j^'Mg. 

We describe the halo density profile by the following form 
lNavarroetaII(ll997h . 



Ph(r|M) 



Ps 



(r/rs)(l+r/rs)2 



(A.8) 



The total halo mass is then written as (see iTakada & JainI 

I2003h : 



M 



-r 

Jo 



4-nrdr ph{r\M) - 



47rps'-J 



ln(l +c)- 



l+c 



, (A.9) 



where c = ryijr^ is the "concentration parameter", for which we 
assume the following expression, 



c(M,z)^ 



Co 

l+z 


M ' 



(A. 10) 



We take cq - 11 and/? — 0.13, and M* is defined such that 
v(z = 0) = 1, i.e. i5c(0) - o-{Mi,). The virial radius rvh is given 
by the following relation 



M = 



47rri 



-poAvii-(z) , 



(A.ll) 
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with Avi,(z) bein g the critical overdensity for virialisa- The two-halo term is derived from the dark-matter power 

tion at redshi f t z ([Kitavama & Sutol Il996t iNakamura & Sutd spectrum and the halo two-point correlation function: 

19971: iHenrvl |2000|). We take the fitting formula from 

Weinberg & Kamionkowsldl (l2003h PiiKr) = Pn,{k) 

dMn(M)-p— 4feh(M, r)\uh(k\M)\ 



Avir(z) = 187:^ 



1+ 0.399 (n„'(z)-l) 



0.941 



(A. 12) 



Jm, 



n' ,{r) 



FoUowing lTinker et alj (l2005h . we use the scale-dependent halo where 
bias 



»;(«..., = .;(«..,"*'"«-" 



1.49 



(A.13) 



«;ai('-) 






M|,„,(r) 



n(M)A^(M,z)dM. 



(A.20) 



(A.21) 



[l+0.69^„(r,z)]209' 

The upper int egration limit Mut^ir) takes into account the 
where ^„^is the matter c orrelation function^ The large-scale halo halo exclusion (!Zhi^'2Q0 |), i.e. the fact that haloes are non- 
overlapping. We follow iTinker et al.l (|2005|) to compute Miini(r) 
by matching n' ^^(r) with the following expression, 



bias bh(M, z) is given bv lSheth et"an (l200ll) as 
1 



bi,{M,z) = ^h(v) = 1 + -^ [ VaCav^) + ^/Eb{av^) 



2^I-e 



y/adc 



(av^y 



iav^y + b(l-c){l-c/2) 



->^ 



(A. 14) 



/ 



As in [Tinker et al.l (|2005|) . we adopt the revised parameters a - 
I/V2, /7 = 0.35 and c = 0.8. 

A.2. The galaxy correlation function 

We write the correlation function as a sum of two components: 
The one-halo term, to express the galaxy correlation inside a 
halo, and the two-halo term, to account for halo-to-halo corre- 
lation. 



dMin(Mi)N(Mi) 

(M2)N{M2)P(r,Mi,M2), 



: fdM2n( 



(A.22) 



where P{r, Mi , M2) is the probability that two ellipsoidal haloes 
of mass Ml and M2, respectively, do not overlap. Defining x = 
r/[rvir(Mi) + rvir(M2)] as the ratio of the halo separation and the 
sumofthe virial radii, and y = (x-0.8)/0.29, iTinker et al.i(i2005l) 
found the probability of non-overlapping haloes to be 



[0 ify<0 

P(r,MuM2) = P(y) = I {3y- - 2y^) if < y < 1 
1 ify>l 



(A.23) 



^(r)=l+^i(r) + ^2W, 



(A.15) 



The one-halo term depends on the number of galaxy pairs 



We Fourier-transform Eq. lA.20l for a range of tabulated values of 
r, to compute the two-halo term ^2 of the galaxy autocorrelation 



<A^(A^ - 1)) per halo. This is comprised of the central-satellite function. Finally, we renormalise it to the total number of galaxy 

contribution (NcNs) and the satellite-satellite term (A^s(A^s - !))■ P^irs: 

Assuming a Poisson distribution, we write r ' r ^ n 



{NAN, - 1))(M) 
(A^cA^sXM) 



N,{M) ; 
N,{M)N,(M) . 



1+^2(0 = 



gnl^ 



"gal 



[l+ftW]. 



(A.24) 



(A. 16) 



The angular two-point correlation function w{0) is computed 
The one-halo term correlation function for central-satellite pairs from the observed photometric redshift distribution and ^(r) us- 
is then ing Limber's equation ( Limber 1954) : 



l+^cs(r,z) = 



dMn{M,z) ; ;^ Ph(r\M) . 



(A.17) w(6') = 2| dxfixf I du^ir^ ^u^ + x^O^), (A.25) 

Jo Jo 



with 



/w = 



n{z) 
dx{z)ldz ' 



The lower integration limit Mvii-(r) is the virial mass contained 
in a halo of radius r, computed with Eq. lA.llI This accounts for 
the fact that less-massive haloes are too small to contribute to the 
correlation at separation r. 

The one-halo satellite-satellite contribution ^ss involves the and x(z), the radial comoving coordinate 
halo profile auto-convolution and is therefore easier to compute 
in Fourier space. The corresponding power spectrum is written 
as 



(A.26) 



PsAk) 
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Jm 



Mhigh 
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gal 



Appendix B: Best-fitting IHOD parameters and 
deduced quantities 



(A.I8) 



where u\i{k\M) is the Fourier transform of the dark-matter halo 
profile ph{r\M). The correlation function ^ss is then obtained via 
a Fourier transform. The one-halo correlation function is the sum 
of the two contributions. 



^lW=^csW+^ssW. 



(A. 19) 
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Table B.l. Description of all galaxy samples and best-fitting HOD parameters. Halo masses are given in h Mq and galaxy number densities in h Mpc 



Redshift 



Mg- 5 log h 



A' log«gal,obs l0gM„ 



log Ml 



log Mo 



0"looM 



log "gain 



logMh, 



;r'/dof 



0.2<z<0.4 <-17.8 451203 

0.2<z<0.4 <-18.8 243 833 

0.2<z<0.4 <-19.8 104 889 

0.2<z<0.4 <-20.8 28 081 

0.2<z<0.4 <-21.8 3 679 

0.4<z<0.6 <-18.8 551547 

0.4<z<0.6 <-19.8 244 854 

0.4<z<0.6 <-20.8 65 557 

0.4<z<0.6 <-21.8 6 724 

0.6 < z < 0.8 < -19.8 458 863 

0.6<z<0.8 <-20.8 141129 

0.6 < z < 0.8 < -21.8 17 683 

0.6 < z < 0.8 < -22.8 656 

0.8<z<1.0 <-20.8 200516 

0.8<z<1.0 <-21.8 31001 

0.8<z< 1.0 <-22.8 1411 

1.0<z< 1.2 <-21.8 52310 
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Table B.2. Description of red galaxy samples and best-fitting HOD parameters. Halo masses are given in h ^Mq and galaxy number densities in h ^ Mpc"* 



Redshift 



Mg- 5 log h 



N log;igai,obs logM„ 



log Ml 



log Mo 



0"logM 



log "gain 



logMh, 



;r'/dof 



0.2<z<0.4 <-17.8 153 729 

0.2<z<0.4 <-18.8 100 829 

0.2<z<0.4 <-19.8 54 393 

0.2 < z < 0.4 < -20.8 17 200 

0.2<z<0.4 <-21.8 1545 

0.4<z<0.6 <-18.8 209469 

0.4<z<0.6 <-19.8 122 073 

0.4 < z < 0.6 < -20.8 42 303 

0.4<z<0.6 <-21.8 5 210 

0.6<z<0.8 <-19.8 216675 

0.6 < z < 0.8 < -20.8 83 687 

0.6 < z < 0.8 < -21.8 11 982 

0.6 < z < 0.8 < -22.8 385 

0.8<z<1.0 <-20.8 111158 

0.8<z<1.0 <-21.8 21494 

0.8<z< 1.0 <-22.8 955 

1.0<z<1.2 <-21.8 30637 
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Appendix C: Two-point correlation function 
measurements 
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Coupon et al.: Galaxy clustering in the CFHTLS Wide 
Table B.3. Description of blue galaxy samples. Galaxy number densities are given in h'^ Mpc-'. 



Redshift Mg — 5 log h N log Ogaiobs 

0.2<z<0.4 <-17.8 297 474 -1.72+°°{ 

0.2<z<0.4 <-18.8 143004 -2.03+;;;°^ 

0.2<z<0.4 <-19.8 50496 -2.49!o;;^ 

0.2<z<0.4 <-20.8 10881 -S-lS^^jJ* 

0.4<z<0.6 <-18.8 342078 -"^-Ql^^of^ 

0.4<z<0.6 <-19.8 122 781 -2.45!;^;j;^ 

0.4 < z < 0.6 < -20.8 23 254 -3.18!;^;;!^ 

0.6 < z < 0.8 < -19.8 242 188 -2.36!0;j;2 

0.6 < z < 0.8 < -20.8 57442 -2.99!0;5;^ 

0.6 < z < 0.8 < -21.8 5 701 -3.99!0« 

0.8 < z < 1.0 < -20.8 89 358 -2.92!;;;;^ 

3 89+°°^ 



0.8<z<1.0 <-21.8 9507 -3.8' 

1.0<z<1.2 <-21.8 21673 -3.62+0 0^ 
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Coupon et al.: Galaxy clustering in the CFHTLS Wide 
Table C.l. Two-point correlation function measurements in the range 0.2 < z < 0.4. 

All galaxies 



e(deg) 


Mj,-51og/i<-17.8 


Mj-51og/i<-18.8 


Mj,-51og/7<-19.8 


Ms-5\ogh<-20.S 


Mj,-51og/i< -21.8 


0.0013 


1.4818 ±0.0322 


2.3199 ±0.0670 


3.5549 ±0.1685 


7.3561 ± 0.9543 


10.3172 ±6.7590 


0.0024 


0.8895 ±0.0213 


1.3596 ±0.0377 


1.9627 ±0.0878 


3.8350 ±0.2778 


7.1094 ±2.8470 


0.0042 


0.5787 ±0.0132 


0.8173 ±0.0222 


1.1466 ±0.0466 


2.0245 ±0.1489 


2.0728 ± 1.2725 


0.0075 


0.3949 ±0.0127 


0.5195 ±0.0180 


0.7166 ±0.0307 


1.0055 ±0.0949 


1.2243 ±0.4861 


0.0133 


0.2819 ±0.0106 


0.3630 ±0.0140 


0.4808 ± 0.0213 


0.7403 ± 0.0500 


0.9581 ±0.2822 


0.0237 


0.1883 ±0.0091 


0.2340 ±0.0130 


0.2893 ± 0.0182 


0.4364 ± 0.0403 


0.5257 ±0.1548 


0.0422 


0.1243 ±0.0072 


0.1501 ±0.0103 


0.1861 ±0.0140 


0.2800 ± 0.0268 


0.4263 ± 0.0780 


0.0750 


0.0776 ± 0.0060 


0.0928 ±0.0081 


0.1 179 ±0.01 11 


0.1700 ±0.0238 


0.1334 ±0.0493 


0.1334 


0.0499 ± 0.0052 


0.0593 ± 0.0068 


0.0761 ± 0.0086 


0.1172 ±0.0173 


0.1055 ±0.0393 


0.2371 


0.0342 ± 0.0052 


0.0426 ± 0.0067 


0.0536 ± 0.0088 


0.0820 ± 0.0169 


0.0888 ± 0.0225 


0.4217 


0.0241 ±0.0046 


0.0297 ± 0.0060 


0.0362 ± 0.0083 


0.0624 ±0.0138 


0.0620 ± 0.0176 


0.7499 


0.0135 ±0.0033 


0.0165 ± 0.0045 


0.0212 ±0.0056 


0.0327 ± 0.0076 


0.0330 ±0.0121 


1.3335 


0.0047 ± 0.0046 


0.0060 ± 0.0047 


0.0082 ± 0.0047 


0.0148 ± 0.0068 


0.0210 ±0.0104 


Red galaxies 


0(deg) 


Mj-51og/i<-17.8 


Mj-51og/i<-18.8 


Mj-51og/7<-19.8 


Mj,-51og/!<-20.8 


Mj-51og/7< -21.8 


0.0013 


4.7974 ±0.1698 


6.0040 ±0.2551 


7.4692 ± 0.5074 


12.0464 ±2.1742 


37.0212 ±31.1297 


0.0024 


3.2359 ±0.0904 


3.7019 ±0.1303 


4.2512 ±0.2044 


7.2855 ± 0.8348 


14.9207 ± 9.6253 


0.0042 


2.1016 ±0.0635 


2.2614 ±0.0780 


2.5039 ±0.1200 


3.5656 ±0.3051 


4.8226 ± 3.5344 


0.0075 


1.4436 ±0.0540 


1.4913 ±0.0594 


1.5948 ± 0.0767 


1.8603 ±0.1843 


4.4201 ± 1.7669 


0.0133 


0.9918 ±0.0456 


1.0089 ±0.0472 


1.0138 ±0.0489 


1.3918 ±0.0982 


3.3844 ±0.9315 


0.0237 


0.6053 ±0.0357 


0.6219 ±0.0381 


0.5959 ± 0.0363 


0.7927 ±0.0671 


2.1519 ±0.4732 


0.0422 


0.3411 ±0.0230 


0.3528 ± 0.0257 


0.3319 ±0.0278 


0.4282 ±0.0418 


1.1750 ±0.2607 


0.0750 


0.1768 ±0.0146 


0.1953 ±0.0173 


0.1961 ±0.0178 


0.2743 ± 0.0344 


0.5895 ±0.1460 


0.1334 


0.0950 ± 0.0098 


0.1094 ±0.0125 


0.1199 ±0.0130 


0.1749 ±0.0205 


0.3728 ±0.1074 


0.2371 


0.0596 ± 0.0092 


0.0728 ±0.0116 


0.0807 ± 0.0127 


0.1151 ±0.0196 


0.2964 ± 0.0609 


0.4217 


0.0409 ± 0.0077 


0.0506 ± 0.0096 


0.0560 ±0.01 15 


0.0847 ± 0.0153 


0.1940 ±0.0458 


0.7499 


0.0239 ± 0.0062 


0.0301 ± 0.0087 


0.0328 ± 0.0090 


0.0462 ±0.01 14 


0.1166 ±0.0292 


1.3335 


0.0101 ±0.0062 


0.0142 ± 0.0086 


0.0142 ± 0.0083 


0.0220 ± 0.01 19 


0.0572 ± 0.0309 


Blue galaxies 


0(deg) 


Mj-51ogA<-17.8 


Mj,-51og/i <-18.8 


Mj-51og/7<-19.8 


Mj,-51og/7<-20.8 


Mj-51og/i< -21.8 


0.0013 


1.0697 ± 0.0363 


1.6188 ±0.0912 


2.4395 ± 0.3425 


5.9775 ±2.1846 


33.1298 ±40.8099 


0.0024 


0.5225 ± 0.0225 


0.7530 ±0.0501 


1.0018 ±0.1375 


1.0833 ±0.5121 


5.9415 ±4.5196 


0.0042 


0.3040 ±0.0145 


0.4330 ± 0.0260 


0.5566 ±0.0585 


1.0238 ±0.2918 


1.3497 ± 1.4446 


0.0075 


0.1801 ±0.0086 


0.2253 ± 0.0130 


0.3712 ±0.0356 


0.6845 ±0.1501 


0.8682 ±0.7125 


0.0133 


0.1339 ±0.0067 


0.1661 ±0.0087 


0.2189 ±0.0229 


0.1885 ±0.0979 


0.2393 ±0.3912 


0.0237 


0.0934 ± 0.0049 


0.1 144 ±0.0078 


0.1490 ±0.0161 


0.1619 ±0.0526 


0.2913 ±0.2202 


0.0422 


0.0681 ±0.0045 


0.0799 ± 0.0062 


0.1063 ±0.0125 


0.1598 ±0.0318 


0.1634 ±0.1230 


0.0750 


0.0507 ± 0.0046 


0.0598 ± 0.0056 


0.0753 ± 0.0105 


0.1043 ±0.0234 


0.0250 ± 0.0562 


0.1334 


0.0369 ± 0.0039 


0.0414 ±0.0050 


0.0560 ± 0.0085 


0.0872 ± 0.0190 


-0.01 17 ±0.0385 


0.2371 


0.0263 ± 0.0039 


0.0301 ± 0.0048 


0.0366 ± 0.0079 


0.0602 ± 0.0165 


0.0452 ±0.0231 


0.4217 


0.0183 ± 0.0035 


0.0214 ±0.0044 


0.0279 ± 0.0079 


0.0443 ± 0.0160 


0.0169 ±0.0182 


0.7499 


0.0098 ± 0.0026 


0.0113 ±0.0030 


0.0178 ± 0.0052 


0.0276 ± 0.0086 


0.0178 ± 0.0104 


1.3335 


0.0029 ± 0.0037 


0.0032 ± 0.0035 


0.0070 ± 0.0038 


0.0133 ±0.0063 


0.01 12 ±0.01 16 
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Table C.2. Two-point correlation function measurements in the range 0.4 < z < 0.6. 

All galaxies 



0(deg) 


M(,-51og/!<-18.8 


M<,-51og/7<-19.8 


Mj-51og/i<-20.8 


M5-51og/i<-21.8 


0.0013 


1.3933 ±0.0262 


2.2367 ± 0.0615 


4.1633 ±0.2645 


9.4751 ±3.9229 


0.0024 


0.8180 ±0.0148 


1.2711 ±0.0286 


2.3764 ±0.1307 


7.9475 ±1.7637 


0.0042 


0.5276 ± 0.0104 


0.7745 ± 0.0193 


1.4069 ±0.0628 


3.8402 ± 0.7920 


0.0075 


0.3481 ±0.0073 


0.4912 ±0.0125 


0.8450 ± 0.0375 


1.8830 ±0.3531 


0.0133 


0.2218 ±0.0055 


0.3019 ± 0.0090 


0.5299 ± 0.0206 


0.8854 ±0.1511 


0.0237 


0.1398 ±0.0047 


0.1862 ±0.0069 


0.3035 ±0.0148 


0.6038 ± 0.0930 


0.0422 


0.0876 ± 0.0035 


0.1167 ±0.0055 


0.1787 ±0.0106 


0.2470 ± 0.0497 


0.0750 


0.0563 ± 0.0029 


0.0749 ± 0.0044 


0.1 150 ±0.0085 


0.2270 ± 0.0288 


0.1334 


0.0378 ± 0.0025 


0.0504 ± 0.0035 


0.0800 ± 0.0052 


0.1409 ±0.0165 


0.2371 


0.0254 ±0.0021 


0.0339 ± 0.0029 


0.0499 ± 0.0046 


0.0948 ±0.0146 


0.4217 


0.0145 ± 0.0017 


0.0201 ± 0.0025 


0.0293 ± 0.0044 


0.0497 ±0.0108 


0.7499 


0.0066 ± 0.0018 


0.0101 ± 0.0027 


0.0153 ± 0.0041 


0.0197 ±0.0081 


1.3335 


0.0022 ± 0.0044 


0.0043 ± 0.0041 


0.0088 ± 0.0040 


0.0124 ±0.0079 


Red galaxies 


0(deg) 


Mj-51og/i<-18.8 


M^-51og/7< -19.8 


Mj-51og/i<-20.8 


Mg-5logh < -21.8 


0.0013 


4.0250 ±0.1 139 


5.0502 ±0.1719 


7.2406 ± 0.5425 


19.5919 ± 9.0027 


0.0024 


2.4772 ± 0.0594 


2.9115 ±0.0828 


4.0217 ±0.2380 


10.5452 ± 2.8842 


0.0042 


1.5544 ±0.0337 


1.7346 ± 0.0467 


2.2900 ±0.1250 


6.5579 ±1.4118 


0.0075 


1.0081 ±0.0254 


1.0944 ±0.0255 


1.3993 ± 0.0555 


3.1910 ±0.4906 


0.0133 


0.6123 ±0.0191 


0.6484 ± 0.0204 


0.8565 ± 0.0352 


1.3759 ±0.2206 


0.0237 


0.3545 ± 0.0142 


0.3719 ±0.0141 


0.4487 ±0.0231 


0.7683 ±0.1319 


0.0422 


0.1899 ±0.0100 


0.2063 ±0.0112 


0.2154 ±0.0189 


0.2399 ± 0.0625 


0.0750 


0.1101 ±0.0063 


0.1223 ±0.0077 


0.1460 ±0.0111 


0.3053 ± 0.0422 


0.1334 


0.0662 ± 0.0049 


0.0749 ± 0.0059 


0.0961 ± 0.0070 


0.1346 ±0.0256 


0.2371 


0.0459 ±0.0051 


0.0500 ± 0.0052 


0.0630 ± 0.0058 


0.1122 ±0.0198 


0.4217 


0.0261 ± 0.0040 


0.0287 ± 0.0042 


0.0366 ± 0.0054 


0.0466 ±0.01 19 


0.7499 


0.0124 ± 0.0037 


0.0139 ± 0.0041 


0.0174 ±0.0051 


0.0242 ± 0.0099 


1.3335 


0.0053 ± 0.0047 


0.0054 ± 0.0049 


0.0075 ± 0.0050 


0.0124 ±0.0090 


Blue galaxies 


0(deg) 


Mj-51og/!<-18.8 


M^-51og/7<-19.8 


Mj-51ogA<-20.8 


Mg-5\ogh < -21.8 


0.0013 


0.9171 ±0.0314 


1.4191 ± 0.0867 


2.4257 ± 0.5698 


-2.0388 ± 2.8547 


0.0024 


0.5111 ±0.0180 


0.6902 ± 0.0447 


0.8761 ±0.2250 


4.0406 ±5.0154 


0.0042 


0.2934 ± 0.0088 


0.4165 ± 0.0244 


0.8823 ±0.1442 


1.9021 ±2.6653 


0.0075 


0.1813 ±0.0068 


0.2587 ± 0.0172 


0.4201 ±0.0649 


-0.0745 ± 0.8562 


0.0133 


0.1165 ±0.0047 


0.1636 ±0.0106 


0.2602 ± 0.0462 


-0.3840 ± 0.4747 


0.0237 


0.0863 ± 0.0027 


0.1086 ±0.0069 


0.1781 ±0.0322 


0.0085 ±0.3271 


0.0422 


0.0610 ±0.0026 


0.0759 ± 0.0043 


0.1388 ±0.0175 


0.2799 ±0.1861 


0.0750 


0.0439 ± 0.0025 


0.0580 ± 0.0034 


0.0942 ± 0.0142 


0.2588 ± 0.0975 


0.1334 


0.0317 ±0.0020 


0.0425 ± 0.0035 


0.0851 ±0.0117 


0.2357 ±0.0651 


0.2371 


0.0206 ± 0.0016 


0.0317 ±0.0026 


0.0593 ± 0.0096 


0.0646 ± 0.0473 


0.4217 


0.0117 ±0.0014 


0.0202 ± 0.0033 


0.0396 ± 0.0094 


0.0974 ± 0.0322 


0.7499 


0.0058 ± 0.0015 


0.0119 ±0.0027 


0.0289 ± 0.0094 


0.0396 ± 0.0289 


1.3335 


0.0025 ± 0.0040 


0.0071 ±0.0039 


0.0219 ±0.0104 


0.0591 ±0.0314 
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Table C.3. Two-point correlation function measurements in the range 0.6 < z < 0.8. 

All galaxies 



0(deg) 


Mj- Slog/; < -19.8 


M<,-51og/7<-20.8 


Mj-51og/i<-21.8 


M5-51og/i<-22.8 


0.0013 


1.5217 ±0.0356 


3.0868 ±0.1218 


8.0221 ± 1.5873 


N/A 


0.0024 


0.8792 ± 0.0158 


1.5273 ± 0.0550 


3.5010 ±0.4573 


N/A 


0.0042 


0.5453 ± 0.0097 


0.9487 ± 0.0349 


1.7392 ±0.2452 


N/A 


0.0075 


0.3274 ± 0.0074 


0.5640 ±0.0190 


0.9377 ±0.1 134 


4.8069 ±4.4791 


0.0133 


0.1943 ±0.0049 


0.3223 ±0.01 13 


0.6846 ± 0.0694 


1.0014 ±1.4080 


0.0237 


0.1227 ±0.0039 


0.1870 ±0.0076 


0.3379 ±0.0336 


1.4934 ±0.8316 


0.0422 


0.0789 ± 0.0036 


0.1129 ±0.0070 


0.1957 ±0.0269 


0.7017 ±0.3584 


0.0750 


0.0537 ± 0.0029 


0.0772 ± 0.0058 


0.1499 ±0.0178 


0.8865 ±0.2413 


0.1334 


0.0375 ± 0.0028 


0.0543 ± 0.0053 


0.0992 ±0.0151 


0.5009 ±0.1653 


0.2371 


0.0252 ± 0.0026 


0.0372 ± 0.0047 


0.0675 ± 0.0129 


0.2492 ± 0.0983 


0.4217 


0.0152 ± 0.0029 


0.0236 ± 0.0048 


0.0478 ±0.0130 


0.2689 ± 0.0748 


0.7499 


0.0070 ± 0.0027 


0.0128 ± 0.0043 


0.0251 ±0.0098 


0.1509 ±0.0683 


1.3335 


0.0037 ± 0.0046 


0.0053 ± 0.0052 


0.0150 ±0.0078 


0.1576 ±0.0655 


Red galaxies 


0(deg) 


Mj-51og/7<-19.8 


M^-51og/7< -20.8 


Mj-51og/i<-21.8 


Mg-5logh< -22.8 


0.0013 


3.6793 ±0.1127 


5.4089 ± 0.3033 


15.8452 ±4.5991 


N/A 


0.0024 


2.1018 ±0.0547 


2.9129 ±0.1258 


6.9176 ±1.0528 


N/A 


0.0042 


1.2737 ±0.0330 


1.7428 ±0.0633 


3.6702 ±0.4636 


N/A 


0.0075 


0.7465 ± 0.0223 


0.9808 ± 0.0354 


1.4866 ±0.2029 


7.6277 ± 8.1243 


0.0133 


0.41 12 ±0.0125 


0.5170 ±0.0179 


1.1126±0.1101 


1.7428 ± 3.0541 


0.0237 


0.2301 ±0.0087 


0.2942 ±0.0116 


0.4482 ± 0.0544 


1.9737 ±1.5615 


0.0422 


0.1334 ±0.0057 


0.1669 ±0.0086 


0.2866 ±0.0401 


-0.1086 ±0.5148 


0.0750 


0.0861 ± 0.0045 


0.1 104 ±0.0071 


0.1700 ±0.0267 


0.5358 ± 0.4698 


0.1334 


0.0590 ± 0.0042 


0.0755 ± 0.0066 


0.1255 ±0.0243 


0.5715 ±0.2197 


0.2371 


0.0407 ± 0.0039 


0.0548 ± 0.0061 


0.0871 ±0.0164 


0.0901 ±0.1194 


0.4217 


0.0244 ± 0.0040 


0.0332 ± 0.0059 


0.0552 ± 0.0146 


0.2193 ±0.1365 


0.7499 


0.0120 ± 0.0042 


0.0175 ± 0.0055 


0.0260 ±0.0101 


0.1040 ±0.0817 


1.3335 


0.0083 ± 0.0061 


0.0096 ± 0.0060 


0.0182 ± 0.0097 


0.0749 ± 0.0386 


Blue galaxies 


0(deg) 


Mj-51og/;<-19.8 


Mg-5\ogh< -20.8 


Mj-51og/i<-21.8 


Mg-5\ogh < -22.8 


0.0013 


0.8487 ± 0.0496 


1.4324 ± 0.2268 


4.1554 ±3.3984 


N/A 


0.0024 


0.4539 ± 0.0264 


0.4523 ±0.1133 


0.8359 ± 0.9450 


N/A 


0.0042 


0.2618 ±0.0144 


0.3612 ±0.0618 


-0.2525 ±0.5016 


N/A 


0.0075 


0.1553 ±0.0083 


0.2681 ±0.0291 


0.0214 ±0.2699 


N/A 


0.0133 


0.1038 ±0.0055 


0.1476 ±0.0189 


0.3495 ±0.1536 


-1.4572 ±0.5388 


0.0237 


0.0771 ±0.0038 


0.1075 ±0.0114 


0.1912 ±0.0780 


-1.1570 ±0.3078 


0.0422 


0.0563 ±0.0031 


0.0834 ± 0.0077 


0.1621 ±0.0520 


1.9834 ±1.5364 


0.0750 


0.0398 ± 0.0024 


0.0624 ± 0.0065 


0.1584 ±0.0294 


1.5887 ±0.7844 


0.1334 


0.0279 ± 0.0023 


0.0431 ±0.0045 


0.1002 ±0.0178 


0.5221 ±0.4101 


0.2371 


0.0194 ±0.0022 


0.0322 ± 0.0042 


0.0503 ± 0.0160 


0.4081 ±0.1886 


0.4217 


0.01 10 ±0.0020 


0.0187 ± 0.0037 


0.0539 ± 0.0144 


0.4225 ±0.1605 


0.7499 


0.0047 ± 0.0015 


0.0094 ± 0.0032 


0.0298 ± 0.0122 


N/A 


1.3335 


0.0010 ± 0.0032 


0.0013 ± 0.0034 


0.0136 ± 0.0078 


N/A 
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Table C.4. Two-point correlation function measurements in the range 0.8 < z < 1 .0. 

All galaxies 



6)(deg) Mg - 5 log h < -20.8 Mg-5 log h < -21.8 Mg-5 log h < -22.8 

0.0013 1.9297 ±0.0672 6.5369 + 0.8511 61.7234 ± 41.9290 

0.0024 1.1410 ±0.0340 2.2368 ± 0.2615 28.1877 ± 19.7251 

0.0042 0.6653 ± 0.0252 1.5711 ± 0.1320 22.8965 ± 9.3422 

0.0075 0.4075 ± 0.0142 0.8834 ± 0.0785 7.9313 ± 3.7334 

0.0133 0.2298 ± 0.0084 0.4598 ± 0.0453 1.5075 ± 1.4777 

0.0237 0.1446 ±0.0068 0.2568 ± 0.0267 0.5253 ± 0.3927 

0.0422 0.1004 ±0.0052 0.1774 ± 0.0162 0.4819 ± 0.1835 

0.0750 0.0715 ±0.0043 0.1192 ± 0.0108 0.5160 ± 0.1432 

0.1334 0.0525 ± 0.0038 0.0933 ± 0.0076 0.2072 ± 0.0657 

0.2371 0.0355 ±0.0041 0.0574 ± 0.0067 0.1539 ± 0.0430 

0.4217 0.0231 ±0.0041 0.0400 ± 0.0066 0.1461 ± 0.0376 

0.7499 0.0156 ±0.0039 0.0269 ± 0.0066 0.1218 ± 0.0368 



1.3335 


0.01 12 ±0.0056 


0.0195 ± 0.0075 


0.1075 ±0.0371 


Red galaxies 


6((deg) 


Mj-51og/i<-20.8 


Mj-51og/i <-21.8 


Mj- Slog/! < -22.8 



0.0013 


4.0563 ±0.1830 


8.0337 ± 1.1816 


62.4523 ±72.7031 


0.0024 


2.1944 ±0.0796 


3.5455 ± 0.4726 


62.2007 ±48.4134 


0.0042 


1.2422 ±0.0465 


2.1323 ±0.2013 


39.5324 ±16.5788 


0.0075 


0.6819 ±0.0246 


1.1544 ±0.1257 


11.4103 ±6.7419 


0.0133 


0.3588 ±0.0178 


0.6279 ± 0.0633 


3.2474 ± 3.0772 


0.0237 


0.2091 ± 0.0127 


0.3088 ±0.0311 


0.4151 ±0.5025 


0.0422 


0.1384 ±0.0083 


0.2302 ± 0.0206 


0.3642 ± 0.2609 


0.0750 


0.0970 ± 0.0058 


0.1410 ±0.0155 


0.5988 ±0.1924 


0.1334 


0.0678 ± 0.0049 


0.1041 ±0.0101 


0.1791 ±0.0805 


0.2371 


0.0460 ± 0.0057 


0.0590 ± 0.0087 


0.1322 ±0.0458 


0.4217 


0.0295 ± 0.0057 


0.0424 ± 0.0068 


0.0712 ±0.0370 


0.7499 


0.0205 ± 0.0057 


0.0267 ± 0.0074 


0.1058 ±0.0316 


1.3335 


0.0157 ± 0.0072 


0.0197 ± 0.0080 


0.0573 ±0.0251 



Blue galaxies 



e(deg) M^-51og/i< -20.8 Mj - 51og/i< -21.8 Mj - 51og/7 < -22.8 

0.0013 1.0353 ± 0.1344 0.8847 ± 1.3392 320.2051 ± 433.1487 

0.0024 0.5535 ± 0.0699 -0.0151 ± 0.7771 123.6391 ± 153.4200 

0.0042 0.2876 ±0.0413 0.6398 ± 0.3317 21.6725 ± 22.7273 

0.0075 0.2044 ±0.0218 0.5090 ± 0.1797 5.3233 ± 5.7498 

0.0133 0.1323 ±0.0117 0.3050 ± 0.1021 1.7491 ± 2.0265 

0.0237 0.0975 ± 0.0086 0.2763 ± 0.0572 2.4051 ± 1.5689 

0.0422 0.0841 ±0.0057 0.1246 ± 0.0291 1.5988 ± 0.7454 

0.0750 0.0543 ± 0.0046 0.0981 ± 0.0227 0.5858 ± 0.3418 

0.1334 0.0430 ± 0.0046 0.0743 ± 0.0126 0.2787 ± 0.1904 

0.2371 0.0304 ± 0.0043 0.0546 ± 0.0106 0.3135 ± 0.1561 

0.4217 0.0209 ± 0.0041 0.0355 ± 0.0068 0.2642 ± 0.0839 

0.7499 0.0135 ±0.0036 0.0262 ± 0.0077 N/A 

1.3335 0.0086 ± 0.0044 0.0182 ± 0.0071 N/A 
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Table C.5. Two-point correlation function measurements in the range 1.0 < z < 1.2. 

M, -Slog /!< -21.8 



0(deg) All galaxies Red galaxies Blue galaxies 

0.0013 3.7142 ±0.2962 6.3867 ± 0.7254 2.2515 + 0.4971 

0.0024 1.6094 ±0.1283 2.9118 ± 0.3126 0.3135 ± 0.3126 

0.0042 0.8737 ± 0.0657 1.7832 ± 0.1386 0.2647 ± 0.1136 

0.0075 0.6176 ±0.0433 1.1093 ± 0.0762 0.2090 ± 0.0738 

0.0133 0.2992 ± 0.0273 0.5005 ± 0.0425 0.1765 ± 0.0437 

0.0237 0.1926 ±0.0158 0.3243 ± 0.0254 0.1098 ± 0.0302 

0.0422 0.1403 ±0.0130 0.2180 ± 0.0208 0.0851 ± 0.0162 

0.0750 0.1003 ±0.0095 0.1545 ± 0.0171 0.0516 ± 0.0137 

0.1334 0.0727 ± 0.0074 0.1027 ± 0.0162 0.0515 ± 0.0089 

0.2371 0.0520 ± 0.0072 0.0742 ± 0.0135 0.0401 ± 0.0081 

0.4217 0.0297 ± 0.0056 0.0428 ± 0.0101 0.0213 ± 0.0068 

0.7499 0.0137 ± 0.0057 0.0190 ± 0.0064 0.0110 ± 0.0068 

1.3335 0.0059 ± 0.0050 0.0076 ± 0.0055 0.0058 ± 0.0054 



Table C.6. Best-fitting parameters (from Eq. IZST i of M^^ and Mi as function of luminosity threshold corrected for passive redshift 
evolution to approximate stellar mass selected samples. Results are given for all and red samples, as function of redshift bins. Halo 
masses are given in /i^'M© 

All galaxies 

Redshift A log(Mi) q-m A log(Mi) q-m 

0.2<z<0.4 0.203 ±0.005 11.354 ±0.012 0.503 ± 0.010 0.274 ±0.017 12.756 ±0.021 0.439 ± 0.035 

0.4<z<0.6 0.301 ±0.195 11.570 ±0.224 0.404 ±0.149 0.269 ±0.181 12.715 ±0.211 0.399 ±0.191 

0.6<z<0.8 0.377 ±0.044 11.717 ±0.063 0.352 ± 0.039 0.514 ±0.002 12.898 ±0.002 0.233 ± 0.002 

0.8<z<1.2 0.345 ±0.006 11.731 ±0.013 0.383 ± 0.006 0.334 ± 0.002 12.759 ±0.005 0.335 ± 0.002 

Red galaxies 

Redshift A log(Mi) om A log(M,) om 

0.2<z<0.4 0.706 ±0.341 12.154 ±0.129 0.307 ±0.161 0.288 ± 0.009 12.930 ±0.011 0.477 ± 0.020 

0.4<z<0.6 0.776 ±0.531 12.209 ±0.198 0.267 ±0.241 0.552 ± 0.007 13.141 ±0.004 0.268 ± 0.008 

0.6<z<0.8 0.761 ±0.553 12.270 ±0.233 0.302 ±0.176 1.033 ±0.114 13.435 ±0.043 0.144 ±0.036 

0.8<z<1.2 0.405 ±0.429 12.037 ±0.502 0.433 ±0.187 0.839 ±0.114 13.426 ±0.083 0.188 ±0.054 
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Fig. C.l. ID (diagonal) and 2D likelihood distributions of best-fitting HOD parameters for the full sample, in the range 0.4 < z < 0.6, 
andMg-51og/!<-19.8. 



34 



Coupon et al.: Galaxy clustering in the CFHTLS Wide 

Table C.7. Best-fitting parameters (from Eqs.l29land[30t of bg and log{Mhaio) as function of redshift. Results are given for all and 
red samples. 

All galaxies 

Galaxy bias b^ Mean halo mass log(Mhaio> 

Redshift abias ^bias Qh h 

0.2 < z < 0.4 

0.4 < z < 0.6 

0.6 < z < 0.8 

0.8<z< 1.2 

Red galaxies 

Galaxy bias b^ Mean halo mass log(Miiaio) 

Redshift abias ^bias Oh h 

0.2 < z < 0.4 

0.4 < z < 0.6 

0.6 < z < 0.8 

0.8<z< 1.2 



.071 ±0.006 


0.211 ±0.011 


13.106 ±0.002 


0.111 ±0.002 


.166 ±0.002 


0.288 ± 0.003 


13.002 ±0.009 


0.185 ±0.009 


.250 ±0.018 


0.327 ± 0.024 


12.846 ±0.011 


0.202 ±0.010 


.430 ± 0.047 


0.391 ±0.030 


12.803 ± 0.001 


0.174 ±0.001 



.331 ±0.011 


0.162 ±0.017 


13.368 ±0.001 


0.073 ± 0.001 


.421 ±0.013 


0.206 ±0.014 


13.246 ±0.002 


0.092 ± 0.002 


.572 ±0.016 


0.243 ±0.021 


13.090 ±0.011 


0.129 ±0.009 


.679 ± 0.049 


0.375 ± 0.027 


12.971 ±0.003 


0.147 ±0.002 
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